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Abstract 

> 

We generalize the Effective Field Theory of Inflation to include additional light scalar degrees of 
freedom that are in their vacuum at the time the modes of interest are crossing the horizon. In 
order to make the scalars light in a natural way we consider the case where they are the Goldstone 
bosons of a global symmetry group or are partially protected by an approximate supersymmetry. 
<^> We write the most general Lagrangian that couples the scalar mode associated to the breaking of 
time translation during inflation to the additional light scalar fields. This Lagrangian is constrained 
by diffeomorphism invariance and the additional symmetries that keep the new scalars light. This 
Lagrangian describes the fluctuations around the time of horizon crossing and it is supplemented with 
a general parameterization describing how the additional fluctuating fields can affect cosmological 
perturbations. We find that multifield inflation can reproduce the non-Gaussianities that can be 
generated in single field inflation but can also give rise to new kinds of non-Gaussianities. We 
find several new three-point function shapes. We show that in multifield inflation it is possible 
to naturally suppress the three-point function making the four-point function the leading source 
of detectable non-Gaussianities. We find that under certain circumstances, ie. if specific shapes 
of non-Gaussianities are detected in the data, one could distinguish between single and multifield 
inflation and sometimes even among the various mechanisms that kept the additional fields light. 



1 Introduction 



Recent cosmological observations are testing the inflationary paradigm in greater and greater 
detail. In this situation it is very useful to describe this epoch using the effective field theory 
approach, i.e. describe it through its lowest dimension operators compatible with the sym- 
metries. This approach in fact allows for a description of a system in the most general terms, 
separating in a clear way what is determined by the UV theory and what is instead just the 
result of the symmetries of the problem. Because this approach concentrates on the theory 
of the fluctuations it is also the approach that most directly connects to what cosmological 
observations are actually testing about inflation. 

The effective field theory of inflation in the case where there is only one relevant degree of 
freedom during the inflationary phase was recently developed in [1, 2, 3, 4, 5, 6, 7, 8]. What 
makes it possible to describe the theory for the fluctuations without any assumption about 
the fundamental degree of freedom that is driving inflation is the fact that the inflationary 
period has to end and give way to the standard FRW cosmology. This simple point implies 
that time-diffeomorphisms are broken and that therefore there is a Goldstone boson associ- 
ated with this symmetry breaking. As typical in the case of Goldstone bosons its Lagrangian 
can be constructed in general terms because it is very constrained by symmetries. In practice 
the Lagrangian for the perturbations in single-clock inflation is constructed by choosing a 
particular time-slicing where the clock field is taken to be uniform. In this frame, the most 
general theory is built with the lowest dimension operators invariant under spatial diffeo- 
morphisms, like g 00 , and K^ u , the extrinsic curvature of constant time surfaces. Invariance 
under time-diffeomorphisms is recovered after reintroducing the Goldstone boson ir, which 
transforms under time diffs. of parameter £°(x, t) non-linearly as ir(x,t) — > ir(x,t) +£°(x,t). 
In the high energy limit it turns out that in most cases one can concentrate on the scalar 
mode 7r which makes the physics appear very transparently. 

This approach to describe inflation has been very successful so far. For example it pro- 
vided a way to explore in generality the possible signatures of single field inflation, identifying 
previously missed ones 1 . Furthermore this approach can be used to translate the constraints 
on the non-Gaussianities obtained from WMAP data directly onto parameters of the La- 
grangian for the fluctuations without any loss of generality [3]. In fact it is the Lagrangian 
for the fluctuations that is directly tested by cosmological observations. This approach of 
using the experimental data to put constraints on the most generic Lagrangian built with 
the lowest dimensional operators compatible with the symmetries is well established in the 
particle physics community, where it goes under the name of Precision Electroweak Tests 
[10, 11], but it is quite new in the cosmological setting. 

The purpose of the present paper is to generalize the construction of the effective field 
theory of single-clock inflation to include additional light degrees of freedom that might be 
present during inflation. This is a complicated task that we will develop in a series of two 
papers. First the Lagrangian for the possible additional light degrees of freedom during 

1 For example, the orthogonal kind of three-point function of the density fluctuations that is currently at 
about 2a level in the WMAP 7year data [9] was identified in this context [3] . 
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inflation will turn out to be much less constrained by the symmetries than the Lagrangian 
for the Goldstone boson of time translations (it). Second, there can be many additional light 
degrees of freedom involved during inflation: scalar, fermions, vectors. Furthermore it is not 
necessary that these additional fields are in the vacuum at the moment of horizon crossing. 
In this first paper we will develop the theory of multifield inflation (or maybe more properly 
multi-degree-of-freedom inflation) for the case where these additional modes are scalar fields 
whose fluctuations are in their vacuum during inflation. We will treat the case of fermions, 
vector fields, and fluctuations not in their adiabatic vacuum in a different paper [12] 2 . 

Even for additional scalar degrees of freedom in their vacuum, the construction of the 
Lagrangian will not be completely straightforward. There are two relevant facts that need to 
be taken into account. First, the scalar degrees of freedom need to be lighter than the Hubble 
scale during inflation in order for them to acquire long-wavelength fluctuations. However it 
is theoretically difficult to have a naturally light scalar: quantum corrections result in large 
contributions to its mass. The only known mechanisms to make a scalar naturally light are 
either supersymmetry or having the scalar fields be the Goldstone bosons (or the pseudo- 
Goldstone bosons) of a global symmetry that is spontaneously broken. We will study both of 
these possibilities in as much generality as we can. It is of course possible that other additional 
symmetries or different representations can be found and their study might lead to interesting 
signatures. What we aim to present here is not just the study of some quite general models, 
but rather a general formalism that can applied to study inflationary fluctuations in multifield 
inflation in complete generality. We further do not address directly the issue of finding UV 
completions to the effective field theories we present here, but we simply stick to technically 
natural low energy effective field theories 3 . 

In the case where the additional scalar degrees of freedom are the Goldstone bosons of 
a spontaneously broken symmetry, following the standard treatment, it is possible to con- 
struct the most generic Lagrangian for the Goldstone bosons with very mild assumptions 
on the symmetry group that is being spontaneously broken. In comparison to the standard 
construction in Minkowski space there will be the important novelty that these Goldstone 
bosons can couple to the Goldstone boson of time translation (i.e. ir). This will allow several 
new terms in the Lagrangian. We will perform the construction both for the case where the 
symmetry group being spontaneously broken is Abelian and also for the case in which it is 
non-Abelian, as there will be relevant differences. We will also discuss effects originating from 
the soft-breaking of these symmetries. 

We will then move on to study the case where the additional light degrees of freedom are 
approximately supersymmetric. During inflation there is a minimum amount of supersymme- 

2 Inflationary fluctuations that are not in their vacuum at the time of horizon crossing have been recently 
considered in [13, 14, 15]. 

3 It is possible that some natural low energy effective field theories do not have a natural UV completion. 
The question of which natural effective field theories admit a natural UV completion is a very interesting, 
general and complicated question. The answer is not known in generality, and we do not address this here. 
We rather stick to the fact that so far in Nature we have found many physical systems that are very well 
described by low energy natural effective field theories. At the end of the day, even General Relativity is 
nothing but a natural effective field theory. 
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try breaking due to the fact that the vacuum energy is non zero. The additional light fields 
interact with the field driving inflation at least gravitationally so the effective energy scale of 
supersymmetry breaking in the sector of the additional fields is at least of order H, where H 
is the Hubble constant during inflation. This is too high an energy scale for the scalars to be 
of interest as they need to be lighter than H to even fluctuate during inflation and a factor of 
at least about fifty below H if one wants the fluctuations they induce to be scale invariant at 
a level consistent with current observational constraints. This would suggest that having an 
approximately supersymmetric Lagrangian is not enough, however the required fine tuning is 
not enormous so perhaps one or a few light fields could remain accidentally. We will study 
this possibility in some detail. 

A second fact that needs to be taken into account in order to describe the role of additional 
light degrees of freedom is that in cosmology we do not directly observe fluctuations of the 
inflaton or the other degrees of freedom. We observe the effect of these fluctuations on the 
cosmological perturbations such as the CMB anisotropies or galaxy clustering statistics. This 
implies that it is important to describe how the fluctuations in the additional light degrees of 
freedom get converted into perturbations of the primordial plasma for example by affecting 
the duration of inflation or by changing the composition of the plasma. In the case of single 
field inflation this relationship is very simple: a fluctuation in 7r amounted to a small time- 
delay in the inflaton trajectory, and therefore in an enhanced expansion of order Hit at the 
end of inflation. These are the so-called curvature perturbations and are the only kind of 
perturbation possible in single-clock inflation because reheating happens in every place in the 
same way leading locally to the same plasma composition. The only difference between the 
various regions of the Universe being how much inflation lasted. 

In multifield inflation the situation is much more complicated. A fluctuation in one of 
the additional scalar fields can be thought of as changing the inflaton trajectory. In order 
to make contact with with cosmological observations one needs to know how much longer 
the new trajectory is compared to the unperturbed one and also if the plasma will have a 
different composition if reheating is reached from a different trajectory. Contrary to the case 
of single field inflation, this piece of information is not fully determined by the Lagrangian 
that describes fluctuations around the time they cross the horizon. How much inflation lasts 
along each trajectory is an integral effect that depends on the entire trajectory and could even 
be determined entirely by events at the end of inflation. The same is true for variations of the 
composition of the plasma after inflation: for example it could be that this is just determined 
at the time of reheating, independently of the Lagrangian for the fluctuations around horizon 
crossing. 

Naively, this fact could stop us from developing an effective field theory for multifiled 
inflation, concluding that a full knowledge of the inflationary Lagrangian in field space is 
necessary in order to determine the density fluctuations. However we will argue that to a good 
approximation all of these effects are generated when the relevant modes are out of the horizon 
and gradients are negligible. This will lead us to conclude that regardless of the mechanism 
of conversion of the additional scalar field fluctuations into cosmological perturbations it 
can be parametrized by a relationship between the perturbations in the plasma and the 
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fluctuations of the additional scalar fields at the time of horizon crossing that is local in 
real space. Fluctuations in our universe are not very large and are approximately Gaussian: 
by Taylor expanding this unknown relationship in small fluctuations we will summarize our 
ignorance about the conversion mechanism with a set of a few unknown numbers. This is 
quite similar to the so-called 5N formalism [16]. Though these parameters can be known only 
with a detailed knowledge of the model, we are able to identify many interesting signatures 
that are independent of these unknown numbers provided that the fluctuations are made 
cosmo logically observable. 

Since the construction of the effective theory for multifield inflation will be complicated we 
start by summarizing the main signatures, many of which are unique to the mulfield inflation 
case. 

Summary of signatures 

• Quadratic Lagrangian. It is possible for the fluctuations associated to the additional 
fields to have speeds of propagation, or speeds of "sound", different from one another 
and different from one. It is also possible to have dispersion relations of the form 
uj ~ k 2 /M, similar to the Ghost condensate case. Just as in the single field case, 
modifying the speed of propagation leads to large interaction terms for the fluctuations. 
If we schematically denote by a one of the additional light scalar fields, then in the 
limit of small speed of sound or non-linear dispersion relation an interaction of form 
<9j7r<9j<T(7 becomes large, although still perturbative. Furthermore in the Abelian case it 
is possible to have a time-kinetic mixing between a and 7r of the form rra. 

• Cubic Lagrangian. At cubic level, in the Abelian case operators of the form & 3 and 
cr((9j(j) 2 produce the same kind of three-point functions that can be present in single- field 
inflation. Furthermore there can be operators induced by the fact that the symmetry 
associated to the Goldstone bosons might have been softly broken at the time. The 
most interesting operator is a{d^a) 2 and it induces a well defined unique shape that 
could be detectable and that is usually associated with a comparable-in-size three-point 
function of the local form and also with a larger-in-size four-point function of local 
form (as we will discuss). Such a complex signature would be a clear indication of 
multifield inflation that is not expected to happen in single-clock inflation. It is also 
possible that operators of the form o~(d<j) 2 , without a Lorentz invariant contraction 
of the indices, might be detectable. In general the presence of the various operators 
associated to the explicit soft symmetry breaking depends on the way that the symmetry 
is explicitly broken. Because of this operators that are naively subleading can become 
more easily detectable or even the leading ones. In the non-Abelian case the situation 
is in general very similar to the Abelian case. However, quite remarkably, for some 
non-Abelian groups symmetries forbid operators of the form & 3 and &(di<j) 2 . In this 
case, the operators coming from the soft-symmetry breaking, if present, are the most 
relevant and a three-point function induced by an operator of the form a{da) 2 could be 
more easily detectable. As in the Abelian case, such a detection would be generically 
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accompanied by a comparable-in-size three-point function of the local form and also 
by a larger-in-size four-point function of local form. If the soft-symmetry breaking is 
small, and the nature of the non-Abelian group is such that the other cubic operators 
are zero, then the four-point function becomes particularly relevant. In the Abelian 
case it is also possible to have a dispersion relation of the form oj oc k 2 with detectable 
non-Gaussianities induced by the operators &(di<7) 2 and (d 2 a)(dia) 2 . This is similar to 
what can happen in single-clock inflation. Both in the Abelian and in the non-Abelian 
cases with oj oc k 2 dispersion relation, we can have a detectable three-point function 
induced by the operators associated to the explicit soft breaking. In addition to all of 
the shapes we have discussed all cases we can have a three-point function of the local 
type generated at reheating. The local kind is the leading three-point function in the 
supersymmetric case. 

• Quartic Lagrangian. At quartic level, in the Abelian case with a linear dispersion rela- 
tion there are three operators & 4 , & 2 (di<j) 2 and (c^cx) 4 that induce a four-point function. 
Notice that only one of these operators, a 4 , is relevant in single-clock inflation [5]. In 
general these operators give an effect smaller than the ones from the cubic Lagrangian 
but we are able identify a set of approximate additional symmetries that can be im- 
posed to make the induced four-point function the leading source of non-Gaussianity. 
These can be either an approximate symmetry a — > —a, or the requirement that the 
Lagrangian is approximately Lorentz- invariant 4 . In this last case we are left with only 
the operator (c^cr) 4 . Notice that this is the first case in which we can have large and 
detectable non-Gaussianities without having a large violation of Lorentz invariance or 
of the continuous shift symmetry in the Lagrangian for the fluctuations. There is also 
another possible local four-point function induced by the operator a 4 associated to the 
small breaking of the Goldstone symmetry. In some region of the parameter space this 
shape can be the leading one but is degenerate with the four point function that can 
be created at reheating. Still, this offers a way of generating a detectable four-point 
function of the local form without at the same time having a comparable or even larger 
three-point function of the local kind. It is also possible to have a dispersion relation 
for the modes of the form oj oc k 2 . In this case the leading interacting operators have 
to come from soft-symmetry breaking, and a local four-point function as generated by 
the a 4 operator, or a four-point function as generated by o 2 {dia) 2 operator m! ight be 
detectable. Notice that in single field inflation with a dispersion relation of the form 
oj oc k 2 there is a larger set of possibilities for operators generating a large four-point 
function [5]. In the non-Abelian case the discussion is very similar to the Abelian case 
with the only difference being that there are additional symmetries that can naturally 
suppress the three-point function making the case for searching for a four-point function 
even stronger. In the supersymmetric case there is another new shape generated by the 
operator a 2 {d^a) 2 which could be marginally detectable and that would be accompanied 

4 In the sense that the scale suppressing operators violating Lorentz invariance is higher that the one 
suppressing Lorentz-invariant ones. 
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by a larger local-four point function. We are tempted to argue that this shape represents 
a smoking gun for supersymmetry as an approximate symmetry during inflation. The 
caveat to this argument is that we find it hard to make statements involving naturalness 
in the supersymmetric case given that an accidental tuning of order fifty is necessary 
as a starting point. Finally, in all of these cases (and in the supersymmetric case even 
without any tuning), it is possible to have a large four-point function of the local kind. 

• Isocurvature fluctuations. In multifield inflation isocurvature fluctuations are possible. 
The results for the adiabatic fluctuations extends quite simply to this case except for 
two differences. One novelty is that now there can be correlations between adiabatic 
and isocurvature fluctuations. In the Abelian case, we can have a three-point function 
induced by operators schematically of the form (da) 3 , and also from others schematically 
of the form a(da) 2 associated to the symmetry breaking. At the level of the four point 
function, the leading one is in general an operator of the form a 4 that induces a local 
four-point function. The story with the non-Abelian and the supersymmetric cases 
proceeds very similarly, with one interesting novelty for the non-Abelian case: it is 
possible to have a detectable four-point function induced by operators schematically of 
the form a 2 (da 2 ) or a(da) 3 whose detection would be a clear indication of the nature of 
the symmetry group being broken. A similar phenomenology holds also in the Abelian 
and in the non-Abelian case when the dispersion relation is of the form u ~ k 2 /M. Of 
course, detectability of these isocurvature effects is strongly dependent on the actual 
size of the isocurvature component which is already constrained to be subleading. 

Distinguishability 

Depending on the details of the situation it is possible that by measuring the departure from 
Gaussianity of the primordial fluctuations we might be able to distinguish between single and 
multifield inflation and even to distinguish among the three mechanisms we have identified 
to keep the additional scalar fluctuations light: an Abelian or non-Abelian global symmetry 
or supersymmetry. This is possible only in the case of detection of some particular signatures 
that Nature might be kind enough to imprint on the cosmological perturbations. 

• There are several observational signatures that would rule out single field inflation in 
favor of multifield inflation. In addition to the observation of isocurvature fluctuations 
and/or a local type of non Gaussianity [2, 17, 18], which have been known as signatures 
of multi-field inflation for some time, there are particular shapes of the three and the 
four-point function that go beyond the local shape and that cannot be generated in 
single field inflation. A detection of a non-Gaussianity induced by an operator we list 
other than those listed in table 2 would be an indication of the multifield nature of 
inflation. 

• If there is a global symmetry, the non-Abelian nature of the group can be inferred by 
detecting a four-point function in the mixed adiabatic-isocurvature fluctuations induced 
by operators of the form a(da) 3 , accompanied by the detection of a three point function 
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of the form (da) 3 in the adiabatic sector, or of the form a 2 (da) 2 without an analogous 
shape in the purely adiabatic fluctuations. 

• Supersymmetry could be identified by detection of a subleading four-point function 
in the adiabatic fluctuations produced by the a 2 (d^a) 2 operator with the dispersion 
relation u ~ k and subleading with respect to a local four-point function by a factor of 
order N e . However in the supersymmetric case the caveats we discussed at the end of 
the paragraph about the quartic Lagrangian in this section apply. 

In Tables 1 and 2 we give a schematic summary of the potential non-Gaussian signals we 
have been able to identify in this paper, together we a summary of the known results from 
single-clock inflation with a continuous shift symmetry. The details are discussed in sec. 5. 

The subject of multifield inflation is very large. Some related earlier papers that focused 
on non-Gaussian 3-point and 4-point functions are for example [19, 20, 21, 22, 23, 24, 25]. 



2 Effective field theory of single-clock inflation 

In this section we briefly review the effective action for single-clock inflation. This effective 
action was developed in [1, 6] and we refer the reader to those papers for a detailed explanation. 
The construction of the effective theory is based on the following consideration. In a quasi de 
Sitter background with only one relevant degree of freedom, there is a privileged spatial slicing 
given by the physical clock which allows us to smoothly connect to a decelerated hot Big Bang 
evolution. The slicing is usually realized by a time evolving scalar <p(t), but this does not need 
necessarily to be the case. To describe perturbations around this solution one can choose a 
gauge where the privileged slicing coincides with surfaces of constant t, i.e. 5<p(x,t) = 0. In 
this 'unitary' gauge there are no explicit scalar perturbations but only metric fluctuations. As 
time diffeomorphisms have been fixed and are not a gauge symmetry anymore, the graviton 
now describes three degrees of freedom: the scalar perturbation has been eaten by the metric. 
One therefore can build the most generic effective action with operators that are functions 
of the metric fluctuations and that are invariant under the linearly-realized time- dependent 
spatial diffeomorphisms. As usual with effective field theories, this can be done in a low energy 
expansion in fluctuations of the fields and derivatives. We obtain the following Lagrangian [1, 
6]: 

Se.h. + s.f. = f d A x ^r g \^-M 2 x R + M^Hg 00 - M^H 2 + H) + 
+ ^M 2 (t) V° + I) 2 + ±M 3 (t) V° + I) 3 + 

-^V + D^, - - ^SK^SK % + ...1 , (1) 



where we denote by 5K^ U the variation of the extrinsic curvature of constant time surfaces 
with respect to the unperturbed FRW: 5K^ U = K^ u — a 2 Hh^ u with h^ u being the induced 
spatial metric, and where M 2j 3 and M 1)2 ,3 represent some time-dependent mass scales. 
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Table 1: Signatures in Multi-field Inflation. In the first column we give the operator generating the 
non-Gaussian signal: operators quartic in the cr's lead to a four-point function, operators cubic in the 
ex's lead to a three-point function. In the second and third columns we explain with which dispersion 
relation the signal can be generated. In the third we explain if the signal can appear in the Adiabatic 
(Ad.) or the Isocurvature (Iso.) fluctuations. In the fourth we state the potential origin of the signal. 
Here Ab. stands for Abelian; non-Ab. stands for non-Abelian, S stands for supersymmetry, and 
R stands for generated by non-linearities at reheating. The subscript s indicates that the term is 
generated by soft-breaking terms. The symbol t represents that such a signal can be generated in the 
case the soft symmetry breaking term is such that it forbids some of the lowest dimensional terms. 
The symbol * represents the fact that the signal is in general subleading, but still possibly detectable. 
In the last column we explicitly mention if the induced signal has a non-vanishing squeezed limit 
and is therefore detectable also in clustering statistics of collapsed objects. 

Let us comment briefly on (1). The first term is the Eistein-Hilbert term. The first 
three terms are the only ones that start linearly in the metric fluctuations. The coefficients 
have been carefully chosen to ensure that when combined the linear terms in the fluctuations 
cancel. The action must start quadratic in the fluctuations. The terms in the second line start 
quadratic in the fluctuations and have no derivatives. The terms in third line represent higher 
derivative terms. Dots represent operators that start at higher order in the perturbations or in 
derivatives. This is the most general action for single field inflation and in fact it is unique [1] . 

The unitary gauge Lagrangian describes three degrees of freedom: the two graviton he- 
licities and a scalar mode. This mode will become explicit after one performs a broken time 
diffeomorphism (Stiickelberg trick) to reintroduce the Goldstone boson which non-linearly 
realizes this symmetry. In analogy with the equivalence theorem for the longitudinal com- 
ponents of a massive gauge boson [26], the physics of the Goldstone decouples from the two 
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Table 2: Signatures in Single-Clock Inflation with a continuous shift symmetry. The signal associ- 
ated to the four-point function when the dispersion relation is of the form uj oc k 2 is quite rich, and 
we refer to [5] for a discussion about it. 



graviton helicities at high enough energies, equivalently the mixing can be neglected. The 
detailed study of [1] shows that in most situations of interest this is indeed the case and one 
can neglect the metric fluctuations 5 . 

As anticipated, we reintroduce the Goldstone boson (ir) by performing a broken time-diff., 
calling the parameter of the transformation — it, and then declaring n to be a field that under 
time diff.s of the form t — » t + £° (x) transforms as 

ir(x) ->■ n(x(x)) = ir(x) -£°(x) . (2) 

In this way diff. invariance is restored at all orders. For example the terms containing g 00 in 
the Lagrangian give rise to the following terms: 

r -+ ^^^ r -+ 9 00 + 29°^+(d.y. (3) 

We refer to [1] for details about this procedure. If we are interested just in effects that are 
not dominated by the mixing with gravity, then we can neglect the metric perturbations and 
just keep the n fluctuations. In this regime, a term of the form g 00 in the unitary gauge 
Lagrangian becomes: 

^oo ^ .^Stt-^ + I^tt) 2 . (4) 

cr 

Further, we can assume that the ix has an approximate continuous shift symmetry, which 
becomes exact in the limit in which the space time is exactly de Sitter [1]. This allows us 
to neglect terms in n without a derivative that are generated by the time dependence of the 
coefficients in (1) 6 . Implementing the above procedure in the Lagrangian of (1), we obtain 
the rather simple result: 



S n = d x 



M 2 l H(n 2 -^) + 2M^(n 2 + ^-n^)-iM^ + ... 



a* I \ a 



, (5) 



5 Equivalently, the neglected effects are suppressed by slow-roll parameters or by powers of H/Mp\. 

6 Notice that this is not always the case. Interesting inflation models both single field and multifield 
have been recently proposed in which the n fluctuations are protected only by an approximate discrete shift 
symmetry. See for example [27, 28, 15, 29, 30]. 



10 



where for simplicity we have neglected the terms originating from the extrinsic curvature as 
they are usually important only in a regime where the space time is very close to de-Sitter 
space [1]. 

We notice that when M 2 is different from zero the speed of sound of the fluctuations is 
different from one. We have the following relationship: 

Mi - • (6) 

There are two independent cubic self-interactions, 7r(<9j7r) 2 and 7r 3 at this order in derivatives, 
which can induce detectable non-Gaussianities in the primordial density perturbations. A 
small speed of sound (i.e. a large M 2 ) forces large self-interactions of the form 7r(<9j7r) 2 , while 
the coefficient of the operator n 3 is not fixed because it also depends on M 3 . Cosmo logical 
data can therefore constrain (or measure) the parameters of the above Lagrangian. This 
approach has been recently applied to the WMAP data in [3], giving constraints on M 2 and 
M 3 , as well as on the higher derivative operators that we have omitted in (5). This is the exact 
analogous of what happens for data from particle accelerators when the precision electroweak 
tests of the Standard Model are carried out [10, 11]. 



3 Additional light scalar fields 

We are now ready to proceed to the construction of the effective action for multifield inflation. 
We consider the case where the only light fields during inflation are the inflaton and some 
scalar fields that we call 07, with 1 = 1, ..,N (as mentioned, we will consider other kind of 
light fields in a subsequent paper [12]). Since quantum corrections typically make it difficult 
to have light scalar fields we consider the only two case we are aware of in which a light 
scalar field can be made technically natural. In the first case the additional fields are the 
Goldstone bosons, or pseudo-Goldstone bosons, associated with the spontaneous breaking of 
a (possibly softly broken) global internal symmetry group. If the broken group is Abelian this 
is equivalent to considering the theory of a scalar with a shift-symmetry. In the second the 
additional light fields are approximately supersymmetric. We start with the case in which 
the additional scalar fields are the Goldstone bosons arising from the spontaneous breaking 
of an Abelian group, then we move to the case of a non- Abelian group, and finally we end up 
this section with the supersymmetric case. 



3.1 Abelian case (TV-shift symmetries) 

We start with the case in which there are N Goldstone bosons arising from the spontaneous 
breaking a U(1) N group. We construct the Lagrangian in unitary gauge. Because of the shift 
symmetry, each 07 appears with a derivative acting on it. This derivative carries an index 
and because of the symmetry under time-dependent spatial diff.s of the Lagrangian in unitary 
gauge, this index has to be contracted in a diff. invariant way or it has to be of the form 
g^d^crj. Thus at leading order in derivatives, there are only two operators to construct the 
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o/'s Lagrangian: g°^d^aj and g^ v d^aid v aj. Keeping in mind that we cannot write tadpole 
terms for the 07's, the resulting Lagrangian has the form: 

Sm.f. = J d 4 x V^g[Mi(t) 2 V° + Wd^i) (7) 

-e^tY^gTd^aj) + e 2 (t) IJ (g ^a I )(g ^a J ) + 
+e 3 (t) IJ (g 00 + W^a^d^j) + e,(t) IJ (g 00 + l)(g^d,ard v aj) + 
+M 2 (*) 2 V° + 1)V^) + 
+Ms(t)-^ IJK (g ^a I )(g°^a J )(g^d fl a K ) + 
+M 4 (t)- 2 ' IJK (g°»d,a I )(g^d,ajd u a K ) + . 

where we have kept only terms up to cubic order in the fluctuations only. Greek indeces go 
from to 3, capitol latin indexes go from 1 to N and the . . . represent terms higher order 
in the fluctuations or with higher number of derivative as well as terms suppressed by the 
small breaking of the shift symmetry of the 07's. e^, with % = 1,...,4, are dimensionless 
time-dependent coefficients, while the Mj's, with i — 1, ... ,4, are time-dependent parameters 
with dimension of mass. It is technically natural to expect the time-dependence of these 
coefficients to be suppressed by the parameters protecting the approximate shift symmetry 
of 7r, the slow-roll parameters. Here we stick to this case although examples where this is not 
satisfied might exist and be interesting. 'The terms in the first line start quadratic in the field 
fluctuations, while the terms in the third and fourth lines start at cubic order. Without loss 
of generality, we can take the kinetic coefficients e[ J to be equal to S IJ , and e 2 to be diagonal 
in field space (e 2 J = e\b IJ \ This can always be done with a proper field redefinition. 

Eq. (7) simplifies when we reintroduce the Goldstone boson of time translations n and we 
go in the high energy regime where we can neglect metric fluctuations 1 . We reintroduce n 
by performing a time diff as in the former section. This time we need the two transformation 
laws: 

a " -> di ^ 1 r -> -m + *) + if^T . (8) 



7 Metric fluctuations can be treated as in single field inflation. In the ADM formalism, one fixes the gauge, 
solves the constraint equations for the constrained variables N, N l , and plugs back the solutions in the action. 
However, the effect from metric perturbations is expected to be very small. If we neglect the term in M(, the 
mixing with gravity starts at non-linear level in the cr's, therefore affecting higher order correlation function. 
This is different to the case of single field inflation where the mixing starts at linear level. Here it becomes 
irrelevant in the limit H/Mp\ very small. The gravitational mixing with the Goldstone boson it is instead 
generically of order of the slow roll parameters. The term in M[ induces a linear mixing, and so it could be 
more relevant, at least in principle. However, as we will see next, in order to avoid ghosts or strong coupling 
in the theory, we need to have M{ < (HM^) 1 ' 4 , which reduces the effect of mixing on non-Gaussianities at 
most of order of the slow roll parameters. 
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We obtain 



5(2,3) 



d A x, 



-9 



Mi 



21 



-2,-^ + ^4 
a 1 



[1 + 71)07 + —di-KdiGi + 



-d^oid^oi + e[ ^-(1 + 77)07 + ^d^dia^j ^-(1 + 77)67 + -^d^dia^j + 



-2e 3 IJ 7T& I & J - 2e 4 7J 7r d^ojd^oj - 4M 2 2J 7r 2 07 + 



-m; 



-2, UK ■ 



-2, UK 



(9) 



where small latin indexes run from 1 to 3. It is useful to join together this Lagrangian and 
the one from single field (5), and to split it into a quadratic and a cubic term. We obtain: 



5 (2) 



(10) 



d A x\/—g 



(2M 2 4 - M 2 x H)k 2 + M 2 X H 



+ 2M{ 1 n& I + (l + e 1 2 )a I & I + 



and 



I 



d 4 xy/^g 



-2Mt tt%£ + ( 2Ml - A -Ml ) tt 3 + 



(11) 



■(Ml + 4M 2 2 ) V6 7 - Ml 1 ^ l7l)2 - - ™ 21 



a 



0/ - 271^ 7T 



2 (e 2 - e 3 + e 4 ) /J 710707 - 2e 4 J 7r 



+ ( M 4 - M 3 1 07070*: - M 4 



-2g 2 - 
07- 



(7/ 



+ . . 



0/ 



In both equations, . . . represent higher derivative terms or terms that break the shift symme- 
try. Let us analyze the quadratic and the cubic Lagrangian separately. 



• Quadratic Lagrangian 

In the 7T Lagrangian the term in (Sg 00 ) 2 induces a speed of sounds different from one for the 7r 
Goldstone boson. Because the Lorentz symmetry is spontaneously broken, a speed of sound 
equal to one is not protected by any symmetry [1]. The same is true for the 07 fields. In 
addition to the standard Lorentz invariant kinetic term for the 07's the operator proportional 
to e2 generates an additional time-kinetic term. This has the effect of changing the speed of 
sound of the 07's: each 07 can have a different speed of sound. This generalizes the case of 
multifield DBI inflation, where it was found that all of the 07 fields and the Goldstone boson 
7T have the same speed of sound [31]. This is a restriction that comes from the particular 
symmetries of the DBI construction but is not general. 

The operator proportional to Mf 1 generates a kinetic mixing between the 07's and ir. 
Again, this mixing is not Lorentz invariant: it is only at the level of the time- kinetic part 
7T07, while the spatial-kinetic part cannot mix through terms of the form diird % ai. This is an 
unexpected constraint that follows from the fact that it non-linearly realizes time diff.s. This 
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is analogous to what happened in single field inflation, where the coefficient of the operator 
7r 2 was not protected by any symmetry, while the coefficient of (diir) 2 was fixed to be equal to 
HMp V Unfortunately this mixing term is not important in a wide range of parameter space. 
In fact, let us consider for simplicity the case in which the mixing is only between n and one 
of the Goldstone bosons 07, that we can call a. In that case upon rescaling a and n as 

M 2 

a res = (1 + e 2 ) V V , 7r rcs = -Mn 71 > ( 12 ) 

and neglecting the time-dependence of the rescaling coefficients, one can put the time-kinetic 
matrix in the form 



<7 r , 

7Tr, 



(0"rcs ^rcs) [-,/-,, x(2M 2 4 -HM2 1 ) J ( . \ = Ores 7T rcs ) ( 

\ 1 (l + e 2 J- ^4 / V TTrcs / V 1 i+eunmi: 

(13) 

where <7 res and 7r rcs represent the rescaled a and 7r fields, and where we have defined e unm i x 
as (1 + e 2 )(2M| — HMp^/Mf — 1. We have also neglected numerical coefficients of order 
one. For e unmix < 0, the time-kinetic matrix has a negative eigenvalue. This means that in 
this limit the theory has a ghost and does not make sense. Clearly in the regime e unm i x 3> 1, 
the effect of the mixing becomes negligible. Therefore, it is only in the region of parameter 
space where < e unmix < 1 that the effect of the mixing is important. As we will underline 
shortly, the kinetic-mixing term is associated with an interaction operator that in the limit 
< e U nmix «C 1 induces a very large level of non-Gaussianity in the density fluctuations. 
So far in the Lagrangian in (7) we have neglected to write down operators of the form 

M~ 2 IJ d 2 G 1 d 2 Gj , (14) 

where two derivatives act on the 07 fields. The reason for this is that generically these 
operators are important only at energy scales of order of the cutoff of the theory and therefore 
negligible at energies of order H. However there is a regime where the operator which contains 
only spatial derivatives can be important. This happens because in inflation we are not 
interested in arbitrary low energy fluctuations but rather in energies of order H [2, 4]. In 
this case (for simplicity we neglect the internal indexes /) taking the limit of very large e 2 
it is possible that for frequencies of order H the spatial-kinetic term of the form (d 2 a) 2 / M 2 
dominates with respect to the standard (fta) 2 term. This happens for e 2 > M 2 /H 2 . In 
this limit the dispertion relation is of the form u 2 ~ k 4 /M 2 , which is the same kind of very 
non-relativistic dispersion relation found in the case of ghost inflation [2, 39]. We see that 
this can also happen even in the case of multifield inflation. 

• Cubic Lagrangian 

The cubic Lagrangian is quite complicated. All the interactions of the form 7r 3 , n 2 a, ira 2 and 
a 3 subject to the constraint that a derivative must act on each fluctuations and the operator 
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must be rotational invariant are present. This leaves us with a number of cubic operators 
equal to 

2 I 3/V I 2 ^ + 1} I N 2 I N{N ±1M± 2 ) | N N ( N + 1) _ 6 + 13iV + 9N 2 + 2iV 3 

2 6 2 3 ' 1 J 

where iV is the number of the a fields (there are 10 operators even for a single additional 
scalar). 

We notice that the coefficients of some of these operators are univocally determined in 
terms of some the operators that appear in the quadratic Lagrangian. This is so because in 
unitary gauge we are allowed to write down operators that are not time diff. invariant. Full 
diff. invariance is recovered once we reintroduce the Goldstone boson ir, but it does so only 
at non-linear level. This means that if at a given order in the fluctuations we have a time-diff 
violating term, upon reinsertion of the n this term will induce operators containing at least 
one extra n fluctuation. 

For example for single field inflaton the coefficient of the operator 7r(<9j7r) 2 was uniquivo- 
cally fixed in terms of the speed of sound and the canonical normalization of the quadratic 
Lagrangian: 

M 2 4 tt 2 ->■ M 2 4 7r(^7r) 2 . (16) 

The same happens here for the non diff. invariant kinetic terms involving the 07 fields. For 
example, the operator e I 2 {g Q ^d IJL ai)' 1 which changes the speed of sound of the fluctuations 
induces the following interaction of the form na 2 : 

(-l + e 2 )a I a I ->■ e 2 — <7j , (17) 

a z 

whose coefficient cannot be altered by any other operator. Notice instead that the coefficient 
of the operator tt& 2 is not fixed by £2 uniquely. 

Similarly, for the time-kinetic mixing 7107, we have 

*}'**, - -MV( ( ^a, + 2i d -^fl] , (18) 

\ gt a z J 

and again the coefficient of these two tt 2 ct operators cannot be changed once M\ is fixed. 
Notice that in the limit e unmix 1 one of the fields becomes very strongly coupled. If e unmix 
is very small, 7r and a are approximately forty-five degrees mixed. Upon diagonalization, one 
of the diagonal fields has a very small coefficient in front of its kinetic term proportional to 
e unmix- After canonical normalization, the operators in (18) represent an interaction for this 
field that scales as 1/e^mix an d the theory becomes strongly coupled in this limit. 

• Quartic Lagrangian 

It is interesting to write down the quartic Lagrangian for the a fields. In unitary gauge, there 
are only three operators which involve four 07 legs that are compatible with the approximate 
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shift symmetry of the 07's. The Lagrangian reads: 

Sl 4) = Jd*x^[M- 4 ' IJKL (g^d,a I d,a J )(g^d,a K d,a L )+ (19) 

fi^.UKL ( g o* ditai ) (g^d.aj) {g^d,a K d,a L ) + 

M~ 4 ' IJKL (g^aj) (g^d.aj) (g°^a K ) (g°»d»a L ) . 

Reinsertion of the 7r will induce some operators that are quintic in the fluctuations, of the 
form 7r<7 4 . We neglect them here. After setting to zero the metric fluctuations, we are left 
with: 

r \ / ~ ~ ~ \ IJKL 

S<$ = I d 4 x^ (Mf 4 - M 2 " 4 + M 3 - 4 J &1&J&K&L+ (20) 

(-2Mf 4 + M^ IJKL &i&j(d l a K d i a L ) + M- 4 ' IJKL (d l a I d i a J )(d i a K d l a L ) 



Notice that the combination proportional to Mi is the only one that is Lorentz invariant. We 
will come back to this point later. We stress that this is not the full Lagrangian at quartic 
order in the fluctuations. We are just considering the terms that are quartic in the 07's. 
In principle, there are also terms like ira 3 , ir 2 a 2 . . . that we have neglected. We did this 
because it will turn out that this quartic Lagrangian in a 4 can generically be important for 
observations inducing a large and detactable four point function. The other terms are less 
important. But in order to explain how this can happen we need to explain how the n and 
the 07 fluctuations are related to observables like the curvature perturbations. We will do 
this in the next section. 

• Soft-breaking Lagrangian 

As we will argue later, the U(1) N symmetry of the Goldstone bosons need not only be spon- 
taneously broken, but also it has to be explicitly broken for otherwise a 07 fluctuation would 
have no effect on the curvature perturbations of the universe. This breaking can be concen- 
trated just around the reheating time, but there might be a soft breaking of the symmetry 
at the time the modes cross the horizon. We will concentrate on this possible breaking in 
this section and come back to the breaking at the time of reheating in the next section. For 
simplicity we will concentrate on a single U(l) symmetry, the generalization to the U(1) N 
case is straightforward. 

Let us imagine that the U(l) symmetry is explicitly broken by a term that transforms 
under the U(l) with charge one, and whose size is controlled by a parameter /i 4 M 4 , with 
M being the typical scale suppressing the higher dimension operators in (7). This implies 
that the soft-symmetry breaking is signaled by the appearance in the Lagrangian of terms like 
e ta which have charge one under the U(l). By treating as a spurion transforming under 
the U(l) with charge minus one, it is possible to construct the potential, which reads: 

y W = - a ^cos(^ + S ) +6l ^cos(2^) + ^ si „(2^) + ... . (21) 
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The terms proportional to &i )2 are higher order in /i/M <C 1, and dots represent even higher 
order terms that are negligible. The effect of the phase 9 is to make a = not a minimum of 
the potential. At leading order in \xjM the minimum is at a ~ M9. Expanding around the 
minimum, with a suitable redefinition of a: 

a = M6 + a , (22) 

we obtain a leading potential given by 

Pleading = fl^/i 4 COS ~ ^/i 4 (l - ^ + 1 ^ + . . .) . (23) 

For the a field to have sizable fluctuations, we need to require its mass to be much smaller 
than H. This implies the condition 

// < H 2 M 2 , (24) 

with H ^ M from imposing the cutoff of the theory be larger than H. Notice that in the 
presence of a mass term to, fluctuations of the a fields keep evolving outside of the horizon 
and decrease their amplitude according to the equation 3H&j + m 2 ai = 0. Since modes that 
exited the horizon at earlier times have longer time to evolve, this effect induces a tilt in the 
two-point function of order m 2 /H 2 . Therefore we need to impose fi 4 /(M 2 H 2 ) < n s - 1< 1. 
The first inequality comes from the fact that the tilt n s — 1 receives additional contributions 
from the scale dependence of H and therefore the actual tilt n s — 1 can be larger than 
fi A /(M 2 H 2 ) 8 . The potential in (21) is even under a — > —a. The leading contributions 
breaking this symmetry come at subleading order in fi/M, and are of the form 

V(a) su b-icading ~ h— sin (2-1 ~ ~ b iH 8 j^ , (25) 

where we have neglected an irrelevant tadpole term that is canceled when we impose the 
minimization at higher order in fi/M. 

The other leading interactions coming from the breaking of the U(l) symmetry are given 
by the operators of the form 

^U(^l-l) - - (26) 



M 2 \ \M J J M 4 M 4 

4 / ~ \ 4 



M 4V ; \M J M 5V ; 

^<da) 2 (cos (I) -l) ^-Ad*) 2 * 2 , 



M 4V ' \ \M ) ) M 6 



8 As we anticipated, we are here interested in the case where the additional scalar fields are directly 
observables and so must have quasi scale invariant fluctuations. We do not study here the case, for example 
as treated in [40], in which the additional fields have mass of order H, do not have scale invariant fluctuations, 
but nevertheless lead to physical effects because back-reacting on tt. We will study these models in a subsequent 
paper. 
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where here (da) 2 stays for (did-) 2 , a 2 or a combination of the two. As we have highlighted, in 
an expanding space the first operator is not a total derivative. Upon integration by parts it 
picks up a factor of H and becomes of the form <r 3 . 

Our assumption of having only one spurion of charge one is technically natural, and it is 
at the basis of the suppression of the a 3 interaction. This will have important observational 
consequences. However, it also possible to add a second spurion /i 4 with charge two and with 
size comparable to //. In this case terms of the form 

P sin ( 2 £) (27) 

are allowed. Upon minimization of the potential, terms cubic in b would survive with coeffi- 
cients of order /i 4 /M 3 : 

" 4siI1 ( 2 s) " P 3 - (28) 

We finally notice that in unitary gauge the parameters /i 4 ,/i 4 , . . . , can depend explicitly 
on time t (for example with a linear dependence). In this case we could even consider terms 
where in unitary gauge we act with a time-derivative on /z 4 , /i 4 , .... For terms linear in /i 4 this 
time derivative could be integrated away by parts but this is not possible to do in general at 
higher order in //, /j 4 , . . . . Upon reinsertion of n the time dependence becomes dependence on 
the combination t + n while every time-derivative in unitary gauge now becomes a space-time 
derivative as for example <9° = g 0p d p . This means that when /z 4 , /j 4 , . . . , are time-dependent, 
we will have additional couplings between 71 and a. Similar considerations apply to the case 
for example when we multiply an operator in a with factors of Sg 00 , the leading being a 
mixing term Sg 00 a. The importance of these additional n — a couplings is quite modest as 
they are already proportional to /i 4 ,/i 4 , .... Further, as we will explain later, they can be of 
some importance only when both n and a are relevant for the curvature perturbations at a 
comparable level. This is in some sense a tuned case whose phenomenology lies beyond the 
scope of the present paper. 

3.2 Non-Abelian case 

In the former section we have considered the simpler case of the Golsdstone bosons 07 's orig- 
inating from the spontaneous breaking of a U(1) N group. We now consider the analogous 
non-Abelian case. While the generalization is straightforward we will obtain some new in- 
teresting signatures that are specific of the non-Abelian nature of the symmetry group. The 
novelty here arises from the fact that if the broken group is non-Abelian the Goldstone bosons 
need not to be protected by a shift symmetry to be light. Therefore we will not be forced to 
have a derivative acting on each of the a fields. At cubic level for example in the Abelian case 
we simply had operators of the form cr 3 and b(dia) 2 . The constraint of having a derivative 
acting on each of the a restricted the possible terms we could write. In the non-Abelian case 
we might have operators of the form a& 2 , a(dia) 2 and a 2 a. In the Abelian case we had these 
operators only after considering terms coming from the soft-symmetry breaking, and their 
form was therefore strongly constrained. In the non-Abelian case these operators can appear 
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already in the limit of exact spontaneously broken symmetry, though the symmetry might 
impose further constraints that we have to explore. 

Let us therefore consider a global symmetry group G that is spontaneously broken to a 
subgroup H. For simplicity, we take G to be compact. The theory for the resulting Goldstone 
bosons together with a historical account of the developments is nicely explained in [44] which 
we follow here. The only generalization we will have to make here is to couple these Goldstone 
bosons to the Goldstone boson of time translations (n). 

Let us set up some notation. We label the generators of H by U, while those of G and not 
of H by x a . The Lie algebra takes the following form 

[ti,t 3 \=iCij k t k (29) 

\Xai •Eb\ ^G a biti -\- lC a bc%c • 

We notice that the broken generators x a 's transform under the same representation of H. 
The Goldstone bosons <r a 's are, apart from a normalization, a parametrization of the right 
coset G/H that we can chose as 

-f({a a (x)}) = Exp [ia a (x)x a ] . (30) 

The Lagrangian for the Gosldstone bosons is the most general Lagrangian built in terms of 
the operator D a ^, which is defined as 

7 _1 (a;)^7(a;) = ix a D afl (x) + itiE ifl (x) . (31) 

At leading order in the fields and in the derivatives, D afl and E iv are given by 

Dan = d » a a + ^CabcOTbd^c + ^ (C cde C bea + C cdi C b ia) v b (T c d^a d + 0(d 3 d^(r) (32) 

X x 

Eifj, = -CabiCTad^ab + -CacdCbdiVaVbd^O-c + O{o O^a) . 

The Lagrangian can also include higher derivative operators which have to be built out of 
covariant derivatives defined as: 

{%D li ) a = d v D a ^ + C iah E iv D hll . (33) 

While the a a transform linearly only under H they do not transform linearly under generic 
G transformation. Instead D afl transforms linearly under all G transformations. In particular 
if we define: 

= D afl x a , (34) 

it transform under an element g G G as 

D'^ = h(a(x),g))D,h(a(x),g)r 1 , (35) 
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where h (a(x),g)) is an element of H that depends on both a and g 9 . 

The Lagrangian for the Goldstone bosons is obtained by constructing the most generic 
Lagrangian with the operators D afl and its covariant derivatives that is invariant under linear 
transformation of the unbroken H group and the remaining symmetries of the problem. This 
automatically ensures invariance under non-linear transformation of the full group G. Notice 
that symmetries dictate that each operator containing Goldstone bosons must contain at 
least one derivative acting on one of the Goldstone bosons. Each term in contains at 
least one derivative, and therefore expanding in corresponds to a low energy expansion. 
The Lagrangian we construct differs from the standard one by the additional couplings to the 
Goldstone boson of time-translations. By now we have learned how to couple to this Goldstone 
boson: we write the action in unitary gauge where we restrict the spacetime symmetries to 
be simply time-dependent spatial diffeomorphisms. The resulting Lagrangian reads 



S a = J d A x^~g Tr [F 2 D^ + F 2 D°D° + F^g 00 + 1)D°+ 

F 2 (g 00 + 1)D^ + F 2 (g 00 + 1)D°D°+ 

FiD^D + F 2 D°D°D + ...] , (38) 

where ^1,2,3,4,5, F 12 are time-dependent numbers with the dimension of mass that we expect 
to be of comparable size. The time-dependence of the coefficients breaks the shift symmetry 
of 7r and is therefore expected to be small, proportional to the slow-roll parameters. We 
have also assumed for simplicity that the generators x a are hermitian (this is true for all 
the Cartan's classical Lie groups) and that the x a transform in an irreducible representation 
of H, the generalization being straightforward. The n gets reinserted in the usual way. In 
particular: 



ft*. (39) 



no _, d(t + n 

dx» 

Upon reinsertion of the field 71, neglecting metric perturbations, we obtain: 

S wa = J d A x^g~ Tr [F^D^D" + F^D°D° + 2F%d^D»D° - 2F^D° + F$(d fl ir) 2 D°+ 

-2Fl-kD^ - 2F^nD°D°+ F^D^D + F 2 D°D°D + . . .] (40) 

where we have kept terms in n only up to those of the form na 2 and n 2 a because we will later 
see that those are the relevant ones. 



9 Using the notation of [44], if g e G we have 

D' ail {x) = ® ab {<j{x),g)D h ^x) . (36) 

where & a b are the matrixes of the adjoint representation of G: 

hxbh- 1 = S> ba x a , (37) 

where h is an element of H, that realize a reducible representation of H. This suggests to conveniently group 
the operators D ap into the matrixes of equation (34) . 
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• Quadratic Lagrangian 



Using (32) we can see that the first two terms in the first line of eq. (40) give rise to the 
kinetic term for the <7 a 's. The first contains the standard Lorentz invariant term while the 
second one contains time-kinetic term. This implies that the speed of sound of the Goldstone 
fluctuations can be different from one: 

c 2 I 2 (41) 

° s (F 2 + F 2 ) ■ W 

As usual, having a speed of sound different from one introduces an interaction with 7r at 
higher order. This comes from the fact that ir non-linearly realizes time-diffs. When the 
speed of sound becomes very small, it is possible that higher- derivative terms, schematically 
of the form 

Tr \9 V D^9 V D^ (42) 

and others with similar index contractions, dominate the spatial part of the dispersion relation 
at horizon crossing, and induce a dispersion relation of the form u 2 oc k 4 as it happens in the 
Abelian case and in the single-clock case. 

The fourth term in the first line denotes a kinetic mixing between o~ a and ir. As in the 
Abelian case, such a term in the quadratic Lagrangian induces interactions between 7r and a. 
In order for this term not to be zero, it needs that Tr[x a ] ^ 0. This in practice signals the 
presence of an Abelian commuting subgroup being spontaneously broken: a kinetic mixing 
between n and a is only possible for Goldstone bosons associated to the breaking of an Abelian 
group. 



• Cubic Lagrangian 

The cubic Lagrangian can be quite constrained in the non-Abelian case, the details depending 
on the particular group being broken. First, the fact that the symmetries do not require a 
derivative acting on each of the fields might allow for interactions of the form a (da) 2 that 
would be the leading ones in a derivative expansion. However, the antisymmetry of the 
structure constants makes these terms vanish. Among the terms with a derivative acting on 
each of the cr's the ones proportional to F\ and F?, give rise to the same kind of interactions 
a(dia) 2 and a 3 that can be produced in the Abelian case. However, the coefficients of these 
interactions are proportional to Tr[x a XfeX c ] which can vanish for some groups (for example 
for SU(2)). Further, there are no interactions of the form rr 2 a in the absence of an Abelian 
commuting subgroup. The only remaining cubic interactions are of the form na 2 , and only a 
small subset of these are present. In fact, while it would not be naively required by symmetries, 
the cr's appearing in the interactions always have a derivative acting on them. In summary, 
in the absence of an Abelian commuting subgroup, we have the following schematic cubic 
interactions: 

&(ditj) 2 , cr 3 , 77077 , -h(di(j) 2 , diTtditra , (43) 
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Just as in the Abelian case the coefficient of diitdi<7 a <7 a is uniquely fixed once the speed of 
sound of the <r's is fixed. Furthermore the first two of these interactions vanish in the case 
Ti[x a XbX c ] = 0. This fact will have important consequences. 

• Quartic Lagrangian 

We first concentrate on the part of the Lagrangian that is purely quartic in the cr's. It is possi- 
ble to write down operators of the form Tt[D°D D°D }, Ti[D^D^D g D\ and Ti^D^D" D v ] 
that give rise to the same quartic interactions we found in the Abelian case. However in some 
cases there are lower-dimensional quartic interactions. The lowest dimensional ones come 
from expanding up to quartic order the terms in the first line in (40). These are schematically 
of the form 

Fi (C ebc C ed a + C d afC b fc + CdaiCbic) X (9 (U (7 a )(7 6 (9 /i (7 c )(7 d , (44) 
F 2 (C ebc C eda + C da fC b f c + C dai C bic ) X ' a O < b 6 " c ' d . 

where we have neglected numerical factors. These quartic interactions are expected to domi- 
nate because they have two less derivatives than the interactions one gets from the operators 
~ which would have a derivative acting on each a. The quartic interactions distinguish 
between Abelian and non-Abelian nature of the broken group. Notice also that in some well 
known groups some of the structure constants are zero and therefore the interactions can be 
further constrained. For example for SU (2) x SU (2) broken to the diagonal subgroup C abc 
is zero. As we will explain in detail later the total antisymmetry of the structure constants 
makes these interactions irrelevant for purely adiabatic fluctuations. The next-to-lowest di- 
mensional quartic interactions come from quartic terms arising from expanding the operators 
proportional to Fi j2 in (40). They are schematically of the form 

Tr[x a x b x c ] C cde da a da b da d a e , (45) 

with two possible combinations of space and time derivatives compatible with rotational 
invariance. It is important to notice that these operators are present only when lower di- 
mensional cubic ones of the form (da) 3 are present, which will make these operators always 
subleading from the observational point of view. We will however see that they can still be 
quite interesting. 

We will find that there exist symmetries that can forbid the leading cubic (and quartic) 
interactions, making interactions of the form ~ D afl 4 the leading ones. These come from the 
Lagrangian 

= J d 4 x^g Tr [dnD^D^D" + ~c 12 D lx D v D»D v + (46) 

5 21 D°D D U D" + 5 22 D°D U D°D U + 
5 31 D°D D°D } 

~ J d 4 xy/^g Tr [x a x b x c x d ] \c xl d il o a d il Obd v o c &' 'a d + c X2 d il o a d v o h d il o c &' 'a d + 

c 2 \0 a 6~bd u o c d l 'a d + c 22 a a d l/ a b a c d u a d + c 3 i& a a b a c &d\ • 
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where £11,12,21,22,31 are time-dependent dimensionless numbers of order unity, and where in the 
second passage we have kept only terms of quartic order in the cr a 's. Notice that, depending 
on the properties of the group, it is possible that the operators multiplying en and £12, as 
well as £ 2 i and £ 2 2, are the same, effectively reducing the number of free parameters. 

• Soft-breaking Lagrangian 

It is quite difficult to discuss in full generality the soft-breaking of a non-Abelian group. 
In order to do this one has to specify the spurionic representation of the terms explicitly 
breaking the symmetry. While in the U(l) case in practice there is only one representation, 
for a non-Abelian group there are many, and the Lagrangian differs in the different cases. It 
is however quite straightforward to realize that all the operators that were associated to the 
explicitly braking in the U(l) for example the a (da) 2 operators, can be obtained also 

in this case. The actual presence (or absence) of each breaking term however depends on the 
spurionic properties of the terms explicitly breaking the symmetry. 

In the case the symmetry is softly broken, the Goldstone bosons have a potential and, de- 
pending on the initial conditions at the beginning of inflation, they can have a time-dependent 
vacuum expectation value (vev). Since a time-dependent vev amounts to spontaneously break- 
ing time-diffs, this means that n will be mixed with those Goldstone bosons a a that acquire 
a vev. Let us take for simplicity the case in which one of the Goldstone bosons a a acquires a 
vev and let us say that this is a"i, as one can always rotate the generators so that this is the 
case. 

In order to construct the Lagrangian for the fluctuations, it is best to proceed in the fol- 
lowing way. We can just take the Lagrangian we have constructed and give a time-dependent 
vev to o\. Since the potential is small we can identify two regimes of interests. The first 
in when the higher-dimensional terms schematically of the form . . . , are assumed to be 
negligible. In this case the equation of motion for the vev, (<7i), will be 

K c > + 3/f K c ) - £ sin (J^j = . (47) 

where the subscript c denotes the canonically normalized field. Here for simplicity we have 
taken all the F's as being the same. Notice that since the inflaton can be given by a com- 
bination of different fields (in principle not even scalar fields), in addition to ci, the Hubble 
constant H is here not univocally determined in terms of /z 4 . Clearly, there are two possible 
sub-regimes for the solution of (47), depending mainly of the value of H. One possibility 
is that the field oscillates around the minimum. In this case the oscillations get damped 
exponentially fast in time. We therefore neglect this possibility. A second regime is when the 
field is slowly rolling down the hill. In this case, we have 

This estimate ensures that the higher dimensional operators D^, . . . , and higher derivative 
operators are negligible. In order to construct the Lagrangian for the fluctuations, we reinsert 
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7r by performing the field redefinition 

<7i(£,t) = (ai)(t + ir(x,t)) + 8au;i(x,t) , (49) 

where 5ai itr represents the scalar field fluctuation remaining after the definition (<7i)(i + ir). 
This happens when the field a 1 is not the only field that develops a time-dependent vev. 
Since the number of light fields is an input in the Lagrangian, the presence of an additional 
clock-field appears as a free choice. By substituting (49) into the Lagrangian and expanding 
in fluctuations, there will be terms independent of the fluctuations or just dependent on ir. 
These should be discarded, as they can just be moved to the part of the Lagrangian that 
has to do with the clock field ir and whose coefficients were generic to start with (in fact, 
there could be additional fields that have a vev and contribute to it). The remaining terms 
involve directly the a a , a ^ 1 fluctuations together with additional couplings involving n 
or 8atr,i- If the 7r field happens to be almost entirely composed of the U\ field, then the 
component ba tT) \ is irrelevant and should be discarded. If instead U\ does not contribute 
relevantly to ir, then the field 5a tr ,i appears in the Lagrangian as the other a a , a ^ 1 
components. Intermediate cases are somewhat tuned and we do not consider them here. 
Similarly if it is dominated by the U\ component, then its Lagrangian will be fully determined 
by the coefficients related to (ai(t)), implying a connection between the coefficients of the 
7r Lagrangian and the Goldstone Lagrangian. However in this case the 7r Lagrangian is not 
very interesting, representing a standard slow roll inflation model. Instead in the opposite 
case in which o\ does not contribute much to ir, then those terms are quite irrelevant for the 
7T Lagrangian. The vev of 0\ leads to additional interactions. The most relevant ones that 
are not present in the absence of symmetry breaking and without the time-depende! nt vev 
a re schematically of the form: 

FTx\x x x a x h \D^DlD h ^ D FTr[x 1 x a x b ](& 1 )C aef (T e d ll (T f d li (T h , (50) 

where the second index in D 1 ' is a spacetime index. Similar combinations of the derivatives 
on the af' b fields are also possible. At quartic level, we have operators of the form 

F Tr^x^D^D^ D F Tr[x lXa x b } (a,) (C cde C fea + C cdi C fia ) ^a c d^ d d^a b , 
Tx^XaX^D^DlD^Dl D Tr^^^C^ r d v a e d v a b ff V , (51) 

and similar index contractions. Notice that the same situation we are describing here in 
the non-Abelian case can happen also in the Abelian case. In that case however the field 
redefinition in (49) is trivial, and we therefore did not treat that case explicitly in the former 
subsection. 

The second regime of interests for the possible background solution of is when the 
speed is very large: ~ F 2 . In this case one should consider the whole set of higher 

dimensional operators, for example a generic function P((o"i) 2 ). In the absence of a large 
gradient in the potential energy as in our case, the only solution that does not redshifts away 
very quickly is that for which P' = 0. If we were dealing with a single scalar field with a shift 
symmetry, this solution would represent the ghost condensate [2, 38, 39] and so we would be 
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tempted to call this case as a 'non-Abelian' ghost condensate. One should consider generic 
functions of Tr[D M D M ], TrlD^D^D^D"], . . . , suppressed by a common scale F with generic 
order one coefficients. However, it seems quite generic, for example by looking already at the 
particular function P(X) = F 4 [a(X/F 2 - l) 2 + (3(X/F 2 - l) 4 + . . .], with X = Tr[D^] 
and a, (3, . . . numerical coefficients, that if &i gets a vev of order F, the additional Goldstone 
bosons receive a mass of order F and can therefore be integrated out of the theory. In this 
case we are back to the standard single-field ghost-condensate which is included in the single 
field Lagrangian 10 . It is conceivable that by choosing a particular symmetry group and/or 
by adding additional symmetries, it might be possible to find cases in which the additional 
Goldstone bosons do not get a mass. The search for such a possibility lies beyond the scope 
of the present paper. 

3.3 Supersymmetric case 

In the former two subsections we have justified the presence of additional naturally-light 
scalar fields by assuming that they were the Goldstone bosons associated to the spontaneous 
braking of some internal symmetry. However it is well known that supersymmetry offers 
another reason for having a naturally light scalar field. In our case the simple fact that the 
universe is inflating spontaneously breaks supersymmetry. In the limit in which we neglect 
the effects coming from supergravity, we can imagine that the sector of the additional fields 
a is decoupled from the inflaton ( i.e. the n field). In this way supersymmetry breaking 
is transmitted to the sector of the a fields at leading order only through the effect that the 
curvature has on the loops [4] . Let us consider for simplicity a simple massless chiral multiplet 
consisting of a scalar doublet a and a Weyl fermion ip a with a superpotential and a Kahler 
potential of the form 

W = AS 3 , K = S f S (52) 

where £ is the chiral multiplet. If we compute loop corrections to the mass of the a field, 
we have that in Minkoswky space the loop from the quartic self-interaction cancels with 
the loop with the fermion contribution so that the mass remains zero. In Minkowsky space 
supersymmetry is not broken. However when we put this same theory in an inflationary 
spacetime and we try to compute loop corrections to the mass we realize that because of 
the inflationary background the scalar and the fermion propagators are altered generically in 
different ways for modes with frequency ou < H . On the other hand for frequencies u 3> H, 
the modes live effectively in Minkowsky space and so they are unchanged [4]. This means 
that loop corrections to the mass will cancel for uj ^> H but not for u < H. This means 
that the natural value of the mass of a is m ~ 5m ~ XH . Even though supersymmetry is 
broken at energies of order H, the effective susy breaking scale for the a field is A H, which 
for A 1 is much smaller than H. This conclusion gets partially altered when we reinsert 
back the supergravity effects. Among the various other corrections, now the scalar potential 
gets multiplied by Exp[i^/M|J. This induces a tree level mass for the a of order H 2 . This is 

10 It is straightforward to generalize this to the case in which multiple degrees of freedom contribute to tt, 
just as did for the case where the vev of & is small (see eq. (48)). 
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bad news for our multifield inflationary model, as we need a mass much smaller than Hubble. 
However, it is still true that loops are cut off at the scale of order H and therefore we can 
imagine to perform the mild tuning to make this tree level mass much smaller than Hubble. In 
practice we require it to be of order (n s — l) 1 ^ 2 !! where n s is the tilt of power spectrum. Loops 
will not renormalize this mass at a relevant level. We conclude that imposing supersymmetry 
in the sector of the a field is a good tool to make the scalar fields naturally light if the sector 
is not directly coupled to the inflaton at the cost of assuming a tuning of order n s — 1 for the 
tree level mass. 

The leading Lagrangian neglecting supergravity corrections reads: 



S a = J [dpad^a* + mW + 3m s A(<7 + a*)aa* + 9\ 2 (aa*) 2 



2 



(a + a*) (cr + cr*) 



(53) 



where we have added a mass term to the superpotential and two of the leading non renor- 
malizable (or irrelevant) terms to the Kahler potential: (a + a*) 3 /M Stl + (a + a* 2 ) 4 /M 2 2 n . 
In the above equation, m s is the supersymmetric mass, of order H 2 in order to cancel the H 2 
mass from supergravity effects. 

If we now consider supergravity corrections, their importance depends on the specific 
inflationary model. For example, in the particular case of small field inflation, the leading 
ones come from the vacuum energy of order H 2 Mp V It is straightforward to check that in this 
case the supergravity corrections to the interaction terms in the superpotential are negligible. 
Instead they are relevant for the interactions associated to the Kahler potential. Symbolically, 
the leading ones are given by 



|Vo*. g <* + *•>*. <«) 

By adding all the possible Kahler terms of the right dimensionality, it is possible to tune away 
these two terms. This would still make the theory technically natural, as loop corrections 
generate those terms at a much smaller value than their naive tree level value. 

Summarizing: in the supersymmetric case we can have the following interactions in addi- 
tion to the ones in (54). In the superpotential we have a cubic and a quartic interaction of 
order: 

\m s a 2 a* + h.c. , A W . (55) 
From the Kahler potential we have derivative interactions of the form: 

^ + a %^a* + h.c. , (gjt^ fygfflV . (56) 



3 n < , * N 4 



M s " ' M 2 



11 Of course here we should include all the interactions allowed by symmetries at the same order of 1/M, 
and they exceeds the two ones we have written. However, as we will see, here we will be interested in just the 
dimensionality of the operators in order to extract the observational consequences of all of them, and the two 
we have considered will be enough. 
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The susy breaking effects generate a cubic term with coefficient of order X 2 H, which is smaller 
than the tree level one coming from the supersymmetric mass m s ~ H. This one is of order 
XH. Notice that, contrary to the Abelian and non-Abelian cases, here interactions without 
derivative acting on any of the a fields are only marginally suppressed by the smallness of 
XH. They are partially protected by supersymmetry. We will see that this will be important 
from the observational point of view. Finally, we stress that all of the above interactions are 
Lorentz invariant. This is so because we have avoided introducing a direct coupling of the a 
to the n Goldstone boson as this would lead to a supersymmetry breaking in this sector at a 
much higher scale. 

Finally, in order to avoid the problem of having to cancel the order H 2 supergravity 
induced mass, it is possible to impose a shift-symmetry on the imaginary part of the chiral 
field a by choosing for example a Kahler of the form 

K = (E + Et) 2 , (57) 

and by forbidding a mass term in the superpotential. In this case the imaginary part of 
the scalar component Im(cx) does not get a mass from supergravity corrections. If the shift 
symmetry of Im(cx) is softly broken the leading interactions involving only Im(o") come from 
the superpotential term. They are of the form A 2 Im(a) 4 . A cubic term proportional to Im(a) 3 
is not generated. Notice that supersymmetry plays an important role in this case. The non- 
renormalization of the superpotential allows to break the shift-symmetry in a more general 
way than in the non-supersymmetric case. 

This concludes the construction of the Effective Lagrangian for the fluctuations around 
the time of horizon crossing. As anticipated, in order to be able to deduce the observational 
consequences of this effective Lagrangian, we need to see how the fields that enter in it are 
related to what we observe in the sky. We are now ready to do this. 

4 Relating to curvature perturbations 

When the quasi de Sitter phase ends, the universe eventually undergoes reheating and reaches 
thermal equilibrium 12 . We will assume for the moment that the final composition of the uni- 
verse is independent of value of the a fields which means that the only difference between the 
various universes is the size of the scale factor. This is what the curvature perturbation ( 
parametrizes. Isocurvature perturbations correspond to the case in which the final composi- 
tion is spatially dependent. We will come back to this possibility later in this subsection. 

In single field inflation, there is a very simple relationship between the Goldstone boson 
7i and the curvature perturbation (. This is due to the fact that in single field inflation 
there is only one classical trajectory, and a n fluctuation corresponds to a time delay on this 
unique trajectory. On the other hand, a ( fluctuation corresponds to how much the universe 
expanded by the end of inflation, on surfaces where the physical clock (or the temperature of 

12 In some models, like the curvaton [20], this might occur after a phase of radiation or matter domination. 
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the universe) is uniform. Because of the uniqueness of the trajectory, it is therefore easy to 
relate n and (. At second order, when the mode is outside the horizon ( is given by [2]: 

( = -Hit + Hir-h + )-Hit 2 (58) 

The linear piece of this relationship is easy to interpret. In the presence of a 7r fluctuation 
that acts as a time-delay 5t ~ n (it has in fact units of time), the universe undergoes an 
enhanced expansion by ( ~ H5t ~ Hn. 



• Abelian case 

In multifield inflation the situation is much more complicated. The novel feature of multifield 
inflation is that there are multiple classical trajectories, associated to the various fluctuations 
of the oi fields. The problem is that it is not easy to relate how a fluctuation in a 07 direction 
occurring about sixty e-foldings before the end of inflation translates into how much the 
universe will have expanded by the end of inflation. At linear level in single field inflation this 
was simply given by Hn because the only possible fluctuation mode was moving along the 
same (unperturbed) trajectory. In general, the relationship between ( and a 07 fluctuation 
will be a complicated function of the entire trajectory, from the time a mode crosses the 
horizon to the reheating time (see Fig. 1). This however will not prevent us from developing 
an effective treatment. 

The question we should answer is how much more the universe expands (or equivalently 
inflation lasts) if we change trajectory by an amount 07. Since we are interested in modes 
that cross the horizon about sixty e-foldings before the end of inflation, the A;-mode associated 
with that fluctuation is outside the horizon for all of those sixty e-folding minus a number of 
order one that corresponds to the time during which the mode crosses the horizon. Let us 
first concentrate on those first order one e-foldings and convince ourselves that no large effect 
is generated during that period. This is so because the 07 fields have an approximate shift 
symmetry during inflation. This suppresses how much the potential can change along the 07 
direction before the end of inflation. 

The only possibilities for a large effect can come from a cumulative effect during all the 
sixty e-foldings or due to a large breaking of the shift symmetry at the time of reheating. In 
fact we expect that at this time all the shift symmetries, both of tt and of the 07 fields can 
be broken. Thus we concentrate on the remaining sixty e-foldings during which the modes 
lie outside of the horizon. In the presence of a 0/ fluctuation larger than the horizon the 
universe in each region of space evolves to a very good approximation as if the <jj field had 
no gradients and it was instead uniform in that region. This means that at linear level no 
matter how complicated the relationship between 07 and ( is it must be a relationship local 
in real space and therefore has to have the form 



COO 



<ti(x) , (59) 





where d(/dai\ is a parameters of dimension of mass -1 which represents the Taylor expansion 
of the relationship between ( and 0/ around the unperturbed trajectory (a I = 0). Notice that 
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Figure 1: Representation of a typical multifield potential. Modes of interest for observation cross 
the horizon about sixty e-foldings before the end of inflation. Therefore, effects coming from the 
evolution of the filds after horizon crossing can be treated locally in real space. Our effective theory 
is more general than this example, as it does not assume that the inflaton is a scalar field. This 
example is however interesting in helping in visualizing the different scales in the problem. 



we have neglected any dependence on 07 at horizon crossing. Note that 07 might not equal 
to zero at that time because of the softly broken shift symmetry. However once a mode 07 is 
out of the horizon it is on an attractor solution: 07 is just a function of 07. We are therefore 
not loosing any information by neglecting terms in 07 in (59). 

This can be extended also to non-linear level. Since from observations we know that the 
non-Gaussianities cannot be very large we can compute them perturbatively as corrections 
to the Gaussian evolution. We concentrate on the three and the four-point function at tree- 
level. In the case of the three point-function two of the fluctuations correspond to observable 
modes and are computed using the Gaussian evolution while the third one is evaluated on the 
first correction to the Gaussian evolution, a term of second order in the Gaussian fluctuations. 
When one takes the expectation value each of the two momenta of the second order corrections 
has to be equal to one of the two momenta of the perturbations which are taken to be Gaussian. 
Since their momenta are out of the horizon this means that all the momenta of interest, even 
the ones that enter in the second order fluctuations are momenta that cross the horizon 
around sixty e-foldings before the end of inflation. Similar conclusions apply for the four 
point function. 

Putting together all of the above reasoning we conclude that even at non-linear level, the 
relationship between a 07 fluctuation and the curvature perturbation ( will be to a good 
approximation a non-linear function that is local in real space and that we can Taylor expand 



da 1 



( u 1 ^ 
a Ax) + 







2! cV/aj 



1_ d 3 ( 
1 J 3! da 1 da jdo~K 



a I (x)a J (x)a K (x) + ... . (60) 
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Here as in (59) ^ 



are the Taylor expansion parameters with dimension of mass n , n 



being the number of 07 fluctuations, and should be evaluated on the unperturbed trajectory 

(*/ = o). 

By including also the contribution from n, we obtain the generic expression 
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tt(x) + 



1 dK 



2! dn 2 
1 



7r(x) 2 + 



a I (x) + 
OK 



(61) 



( \ ( \M 1 



3! doidojdoK 



2! <97r<90 7 
cr/(a;)crj(x)crA:(x) + ... . 



<Ji{x)<Jj(x) + 



The terms in the first line represent the linear relationship between ( and n and 07, and they 
are important to give the overall normalization of the power spectrum. The terms in the 
second line are instead important for the three-point function, while the terms in the fourth 
line (of which we included only the term that we find can be more important observationally) 
are important in the four-point function. 

Notice that the terms in the above Taylor expansion that contain only ir can be read off 
from eq. (58): their overall mass scale is given by Hubble. They are fixed by symmetries and 
are not model dependent. The coefficients proportional to 07 are instead model dependent, 
and their magnitude represents how much a 07 fluctuation affects the overall duration of 
inflation. When these coefficients are large enough so that the 07 fluctuations are as or more 
important for ( than the n fluctuations, then we are in a truly multifield inflationary model. 

In the second line the term mixed in ir and aj deserves a special discussion. It represents 
how much a 07 fluctuation affects the subsequent 7r fluctuation (or viceversa). Clearly, because 
of the approximate shift symmetry around horizon crossing of the 07 fields, this effect is 
suppressed by the breaking of the shift symmetry either of the n field or of the 0/ field, and 
can therefore be neglected. 

In eq. (61), the field fluctuations should be thought of as evaluated at the time when the 
modes of interest have just crossed the horizon. Neither the 7r fluctuations nor the 07's are 
constant outside of the horizon. However one can show that for crj = 0, the combination 
in (58) is constant at non linear level when gradients are negligible [2, 17]. The situation is 
a bit more complicated for the 0/ fluctuations. They are not constant, and they get some 
non-linear corrections from the evolution outside of the horizon. We can deal with this by 
treating the fluctuations of a I after horizon crossing as initial conditions that will determine 
how much the universe will expand given such an initial condition. In addition because of 
the non-linear evolution outside of the horizon, the initial distribution of the 0/ fields will 
receive a non-Guassian correction (i.e. even an initially Gaussian distribution becomes non- 
Gaussian). However as noted in [32] since this evolution will occur outside of the horizon and 
can be treated perturbatively the final form of the correction will be local in real space: 



■(x) = a T (x) + /nl'coI; 1 Mx)a K (x) - (<jj(x)a K (x)) + ■■■) 



(62) 
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Naively, this correction is expected to be small because even though it comes from the build 
up during the sixty e-foldings of time of evolution outside of the horizon, the non-derivative 
self couplings of the 07 fields are suppressed by the parameter controlling the breaking of the 
shift symmetry for most of the duration of inflation and therefore the final effect is suppressed. 
However in the last part of the evolution the shift symmetry of the 07 and n fields may be 
violently broken which could make the effect more relevant. 

Because of the local structure of the induced non-Gaussian correction, implementing 
eq. (62) amounts just to a redefinition of the coefficients that involve derivatives with re- 
spect to 07. Since we do not know them the result in (61) is unchanged. Notice that this 
redefinition does not effect the overall suppression of the term in naj in (61). 

There is one final subtlety to notice about the former treatment. Once we look at the 
7r and 07 fluctuations in Fourier space, the various modes will have crossed the horizon at 
slightly different times (see Fig. 1). This means that we have to allow the various Taylor 
coefficients as well as the amplitude of the modes to have a slight dependence on the time 
of horizon crossing. This will allow us to recover information about deviations from scale 
invariance. This is very similar to what happens for single field inflation [1], and the situation 
is treated analogously. The only difference is that while in single field inflation the tilt was 
just a function of the parameters of the unperturbed trajectory (H, H, and H) and of the 
speed of sound of the 7r fluctuations here the tilt is allowed to be a more general function. 

Our relationship (61), even after the redefinition (62), is quite similar to what is usually 
referred to as the 5N formalism [16]. However, the redefinition in (62) shows that our pa- 
rameterization, at least for the quadratic terms, is different: for example, the coefficient of 
the 070 j term is not equal to the second derivative of the overall expansion with respect to 
the 07 fields. However this difference only enters at non-linear level. Furthermore we stress 
again that this approach neglects the contribution to the expansion coming from the order 
one e-foldings during which the modes cross the horizon, and that we argue is negligible at 
leading order. The literature making use of the SN formalism is very vast. For example it 
has been applied in many models of multifield inflation to calculate the cumulative effect on 
the two-point function of the different inflationary trajectory (see for example [16]). After a 
generalization to treat non-linear effects the formalism was applied to higher order statistics 
(see for example [24, 21, 25]). There are also models where the effects on the overall expan- 
sion originates mainly from variation in the reheating process [32] or from a new late phase 
of matter domination as in the curvaton model [20]. 

• Non-Abelian and Supersymmetric cases 

The relationship between ( and the Goldstone's boson 07's in the non-Abelian case proceed 
as in the Abelian case except that in the non-Abelian case the Goldstone bosons transform 
under some representation of the unbroken group H . Let us consider the linear term in <7j 
from eq. (61). First of all in order for there to be a linear term the group H needs to have 
been broken by the time of rehating. This is a reasonable assumption, as the full group G 
might have been softly broken already at the beginning of inflation. Let us consider the case 
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in which the Goldstone bosons transform in an irreducible representation of the unbroken 
group H as we did in (38). In this case, we could perform the following transformation. 
Given h G H, we perform 



da! 



cri{x) 



da K 

where we have chosen the particular transformation such that 



(63) 



da K 

and we have redefined 



6n , (64) 
o 



a\ = 2}{h) LJ aj{x) . (65) 

Independently of the fact that the group H might be strongly broken by the end of inflation 
the transformation o — > a' leaves invariant the Lagrangian in (38), with the exception of 
the soft-breaking that might be present at the time when the modes cross the horizon which 
we keep general in any case. This means that at the level of the linear dependence of ( on 
the <7i fields we can take ( to depend on a single representative mode for each irreducible 
representation of the group H. The quadratic and higher order terms in the relationship 
(61) are unaffected by the transformation in (63) and are therefore generic. This will be 
particularly relevant for the generation of non-Gaussianities. 

The case where the additional light fields are approximately supersymmetric proceeds in 
exactly the same way with the only different that now ( has to depend both on the real and 
the imaginary part of the a fields. Notice that the fact that the a fields have an influence 
on the curvature perturbations means that there is a coupling with the inflaton that breaks 
supersymmetry directly at some lower energy also in the a sector. 

In summary: curvature perturbations are given by (61) where the fields are evaluated at 
a time just after the observable modes have crossed the horizon during inflation. While the 
coefficients of the terms in 7r are at most of order H, the coefficients of the terms involving only 
powers of the cr/s are instead free. These coefficients effectively parametrize our ignorance 
about the full evolution of the fields during the sixty e-foldings of inflation and the epoch of 
reheating. As we let them vary, we interpolate from effectively single field inflation to a truly 
multifield inflationary model. Contrary to the case of single field inflation, this treatment 
cannot be made indefinitely exact. The error is associated with the fact that eq. (61) does 
not properly account for those order one e-foldings during which the 07 fields are crossing the 
horizon. There is an irreducible error associated with this effect that we expect to be of order 
1/N e , where N e ~ 60 is the number of e-foldings that inflation lasts since the scales we observe 
left the horizon. In reality we expect the error to be even smaller because of the 07 fields are 
approximately massless and that should make only the last e-foldings of evolution relevant 
for (61). Since the term in nai is expected to be small, and since as we will explain later 
we expect the term in <ji<jj<jk also to be irrelevant, we are left with a number of unknown 
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coefficients equal to 



N r (N r + 1) 
2 




non — Abelian case , 



Abelian case , 



(66) 



r 



r 



2N + N(N + l) = N(N + 3) , 



supersymmetric case , 



where for the non- Abelian case the sum of r runs over the number of reducible representations 
of H each one with dimension N r , and for the supersymmetric case we have assumed that 
there are N chiral multiplets unrelated by any symmetry. 

It is worth to explicitly point out that in order for the additional fields to play a relevant 
role for the cosmological perturbations, their internal symmetry (either the global one or the 
supersymmetric one) needs to be explicitly broken. As we have seen, there are small breaking 
effects already in the Lagrangian that describes the fluctuations around horizon crossing, but 
there can be additional effects that become relevant around the time of reheating or even later. 
The actual way the symmetry is broken at this stage only has one more relevant observable 
consequence, determining explicitly the coefficients of the terms in (61). Determining their 
size would require a complete description of the model both during inflation and during 
reheating and it probably faces, in each specific model, additional constraints coming from 
the explicit breaking of the internal symmetries in each model. Such a description seems to us 
difficult to achieve in more specific and yet still general terms than the ones we have provided 
here. It lies beyond the scope of the present paper. 

4.1 Isocurvature fluctuations 

It is possible that the composition of the plasma after reheating depends on the particular 
inflaton trajectory, i.e. on the value of the o/'s. Fluctuations of the composition of the plasma 
on surfaces of constant value of the physical clock are known as isocurvature fluctuations. 

Given the discussion of the former subsection, it is straightforward to offer an effective 
parameterization for the isocurvature fluctuations. The composition of the plasma can be 
parameterized by the various ratios of the densities 13 . Let us denote by n« the number 
density of some species i 



where the index i runs over baryonic, neutrino and all the possibile other species. The number 
density of photons is denoted by n 7 . The composition (67) is expected to be a function of 
some coupling constant and decay channels that could depend on the inflaton trajectory. 
Since the interesting modes are all well outside the horizon these functions will be local in 

13 There are other isocurvature perturbations related to velocities and higher order moments of the distribu- 
tion function of the various components of the primordial plasma [45] . They can be included in our treatment 
in a straightforward way but for clarity in the main text we focus explicitly on perturbations in the plasma 
composition. 
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real space and thus can be Taylor expanded around a I = as we did for (61). The additional 
information we need is that the perturbations to the composition of the plasma have to go to 
zero in the limit in which the 07's do: a simple n fluctuation is just a time delay of the same 
unperturbed solution and therefore reheating in the presence of only a 71 fluctuation happens 
in the same way as in the unperturbed case. We obtain: 



1 9(rij/n 7 ) 

a I (x) + 



2 da id a j 



oi{x)oj{x) + 



dndaf 



7r(a;)o"/(a;) . 

(68) 



The same arguments we used in the former section imply that the coefficients of the terms in 
itai are suppressed by the breaking of the shift symmetries and is therefore small. 

In the non-Abelian case we can perform a rotation such that ( depends only on one 
particular 07. However this need not be the same one as the one for the adiabatic fluctuations. 
This will be important for the non-Gaussianities. 

We are left with a number of model dependent isocurvature parameters equal to 

nx (^N + N ( N ^ + = ^N(N + 3) , Abelian case , (69) 

J]?ix(l + N r( N r + !) ) = ^ ^(jV r 2 + N r + 2) , non - Abelian case , (70) 



2 ' ^ 2 

r 

n x (2JV + N(N + 1)) N(N + 3) , supersymmetric case 



(71) 



where n is the number of isocurvature fluctuations we wish to consider. These are the pa- 
rameters constrained by experiments. 



5 Signatures 

5.1 Detectable four-point function from multi-field inflation 

We now begin the analysis of the signatures of multifield inflation in light of the effective field 
theory we have developed. We start with the four-point function. In Ref. [5] it was shown 
that in single field inflation it is possible to have a large and detectable four-point function 
without having at the same time a large and detectable three-point function. This was due 
to a parity symmetry n — > —n that turned out to be an accidental approximate symmetry 
once we introduced a large coefficient for the operator (Sg 00 ) 4 . This had the consequence of 
allowing a large coefficient for the operator 7r 4 without at the same time a large coefficient 
for the (dn) 3 operators. One can therefore have a large four-point function without a large 
three point function. In the case of a linear dispersion relation u oc k only one shape for the 
four-point function is possible, while in the case of a non-linear dispersion relation u oc k 2 
many possible shapes are allowed [5] (see Table 2). 
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In this section we are going to show that in multifield inflation it is possible to have a 
larger class of four-point functions that can be observationally relevant. In fact contrary 
to the case of single field inflation it is possible to impose symmetries that shut down all 
the cubic self-interactions and also those quartic in 7r without shutting down the cr's quartic 
self-interactions 14 . This makes the four-point function the leading source of non-Gaussianity 
in these models. In the case there is a large cubic self-interaction, there is a non-vanishing 
four-point function. However, the signal in the three-point function is much larger [5]. 

There are at least two symmetries that shut down the interactions that involve a n field. 
One is of course to impose a 7r-parity n — > —n to be a good approximate symmetry. This 
shuts down all the cubic interactions, including the ones in na 2 . Notice that the operator 7r 4 
that was studied in Ref. [5] becomes more and more suppressed in the limit in which we make 
the 7r — > — 7r symmetry more exact, as it is originated by the operator (Sg 00 ) 4 which does 
not respect this symmetry. However the parity symmetry for n can never be exact. There 
is a minimum violation originating from the time dependence of the Hubble constant. Upon 
reinsertion of ix this induces terms with all powers in it. This implies that the breaking is 
suppressed by slow roll parameters which in turn means that the resulting cubic operators 
induce a negligible three-point function. For the U(1) N case, by adding an approximate 
oi — > —oj parity symmetry (that for example can become strongly broken only at the time 
of reheating), all the cubic interactions are suppressed. In the non-Abelian case the a\ are 
automatically suppressed. 

Another way to shut down all the cubic interactions is by imposing that Lorentz invariance 
be a good approximate symmetry for the fluctuations. In fact at least in single field inflation 
with an approximate continuous shift symmetry for the Goldstone boson, the only way one 
can have a large three-point function is for the 7r Lagrangian to be very non-Lorentz invariant. 
For example, we saw that large non-Gaussianities were obtained in the limit of a very small 
speed of sound or in the case the 7r 3 term was made large. In multifield inflation we can 
instead require that the scale suppressing the Lorentz- violating operators be much higher than 
the one suppressing the Lorentz-invariant ones. This effectively makes Lorentz symmetry an 
approximate symmetry for the flucutations. This symmetry removes all the n self-interactions, 
and it also removes all the remaining cubic interactions in cr 3 , in n 2 a and na 2 . 

Neither of the symmetries just mentioned shuts down interactions of the form (<9cr) 4 which 
can lead to large and detectable four-point functions. Before proceeding to study them we 
point out that in order for the non-Gaussianity induced by the cr's to be relevant, the cr's 
need to play a relevant role in determining an observable quantity. For example, this is the 
case if the influence of 7T on ( is much smaller than the one from the cr's: 



where cj is the speed of sound of ( and cr fluctuations respectively. Here for simplicity we 
have considered the case in which all the fields have dispersion relations of the form oj ~ c s k, 

14 In this section we will neglect to specify to indexes of the a fields expect for the cases where it is required 
for clarity. 
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a generalization to the case where some of the fields have a dispersion relation of the form 
lu ~ k 2 /M is straightforward (see [2, 3] for the size of the fluctuations in this case). 

There is at this point a very straightforward way in which the a fields could induce a large 
( four-point function, while inducing a negligible three-point function: this is by imagining 
that the coefficient d 3 (/da 3 in (61) is much more important than the coefficient d 2 (/da 2 . 
The term in d 3 (/da 3 induces a four-point function that we call local and we paremetrized by 
g^L- The name local comes from the fact that the four-point function in this case is induced 
by a local-in-space non-linear relationship between ( and 07. Its size is of the order of 
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Analogously, d 2 (/da 2 induces a local three point function parametrized by f 1 ^ of the form 
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In order for the four-point function induced in this way to be observationally larger than the 
three point function, we need that 



9nl( ^ }nl 



d 3 ( 



8a 3 



> 



1 d\ 



C da 2 



da 



(75) 



Now, the derivatives d n (/da n represent how much the expansion of the universe depends on 
the various trajectories, parametrized by the values of o. Their presence is associated with 
the breaking of the shift symmetry along the a directions that can happen around the time 
of reheating or as a cumulative effect during the approximately sixty e-foldings of inflation. 
The important point here is that once the linear term is generated with a scale of order 



da 



1 



(76) 



where M c is a mass scale associated with the conversion mechanism, we expect all the higher 
derivatives to scale as 



1 1 d n -\ 



a M r da 



n-l 



(77) 



where a is a numerical coefficient representing the efficiency of the conversion mechanism. 
The dependence on a can arise when ( is a local function of a given by the following form 



C = OL Ci 



a 



+ c 2 



a 



+ c 3 



a 



+ ... 



(78) 



where ci i2 ,3 are order one coefficients. M b represents a scale associated to the breaking of the 
symmetry protecting the lightness of the a fields at the scale of reheating, and we expect 
all powers of a to be suppressed by the same scale M b . The overall conversion of the a 
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fluctuations into ( can be controlled by an additional parameter a, that we can imagine as 
representing the inefficiency of the conversion mechanism, for a < 1, or the efficiency in the 
opposite case (see [32] for a first realization of this possibility). Given the current constraint 
on the local J'nl, we have that a is larger than order lO^-lO -2 . This means that for the 
relationship in (75) to be satisfied, a tuning at least of about a/( ~ 10 3 is necessary. That is 
to say we do not expect the four-point function generated this way at the time of reheating 
to be more relevant than the three-point one. This leads us to consider other ways in which 
a four-point function can be made larger than a three-point function. Let us start with the 
Abelian case. 



• Abelian case 

One way to make the quartic interactions the leading ones is by imposing that Lorentz invari- 
ance be an approximate symmetry of the Lagrangian for the fluctuations, in the sense that 
the scale suppressing the operators that violate Lorentz invariance is made parametrically 
higher that the scale suppressing the operators that respect it. This makes irrelevant all the 
7r self-interactions, the mixing term in 7ia, all the cubic interactions involving any power of 
the a's (with one derivative acting on each one) and the 7r's, and of all the quartic interactions 
in (20) but the unique combination that is Lorentz invariant: 

M%ad^ad u ad u a , (79) 

where for simplicity we are suppressing all the internal indexes. In order for the four-point 
function produced by this operator to be relevant, it is not enough for the 07's to be non- 
Gaussian, they also need to play a relevant role in determining an observable quantity. For 
example, we could ask that they dominate £, which requires eq. (72) to be satisfied (in this 
case with c s — 1). The interaction in (79) induces a qnl of the form 



A 2 A 
9nlQ ~ -y- 



i? 2 <7 2 1 t o m H 4 

~— g NL ~W 10 —, (80) 

e~h M 4 M 4 fj| M 4 



where we have used that 



c ~ 10- 5 (K 



daj 



H . (81) 





The resulting value of q^l can be very large. The cutoff of the theory is of order M, and 
therefore qnl is just bounded to be smaller than 10 10 . In particular £4/ ' C2\e~h can De l ar g er 
than 10" 4 and thus detectable. In determining the thresholds for detection, we will estimate 
the size of the constraints as [34] 



MInlC] ~ , A[^ L C 2 ] ~ , (82) 

pix pix 



where N p - lx is the number of signal dominated modes. Given the current experimental bounds 
we will impose the thresholds for detection by imposing the above ratios to be tentatively 
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larger than ~ 1CT 4 . Given the fact that we will neglect all the numerical factors, we will 
interpret these limits loosely. We leave a detailed study of the quantitative predictions and 
of the detectability thresholds, along the lines of what done in [3] for single-clock inflation, to 
future work. 

At the level of the cubic operators, the requirement of having quasi-Lorentz invariance has 
got rid of the cubic operators with one derivative acting on each fluctuation. This symmetry 
does not forbid operators with one more derivative, such as Da(d fJi a) 2 /M 3 . If these operators 
were to contribute to the three-point function, their signal would be larger than the one on the 
four-point function. However, it is quite straightforward to realize that operators involving 
□cr do not contribute to the three-point function at tree-level, as they are proportional to the 
linear equations of motion, by which Oa = 0. Operators of the form V '^d^ad 1 'ad^a can be put 
in the same form by an integration by parts. In an expansion in derivatives, the next cubic 
operators allowed by this symmetry contain six derivatives and are therefore subleading. 

We would like to stress that imposing the Lagrangian to be quasi Lorentz invariant has 
reduced the large number of interaction operators to just one. There is therefore only one 
shape for the four-point function in this limit, which is a smoking gun for multifield inflation 
and strong evidence for approximate Lorentz symmetry of the Lagrangian of the inflationary 
perturbations 15 . 

Notice that the quartic operator in (23), schematically of the form cr 4 , that is present in 
the case of explicit symmetry breaking, is not eliminated by imposing an approximate Lorentz 
invariant, and his contribution can be important. Comparing with the leading term we are 
considering we obtain: 



Here we have normalized the cr terms in the symmetry breaking potential with a scale A^g 
that represents the unitarity bound coming from the interactions associated to the symmetry 
breaking potential. The scales A^g and M are independent (M is the the unitarity bound 
coming from interactions compatible with the shift-symmetry). Because of the discrete shift- 
symmetry cr — >■ a + 27rAu,s, the interactions suppressed by M do not renormalize A^g and 
because of the continuous shift symmetry they do not renormalize the scale [i either. On the 
other hand, interactions in /i 4 generate derivative interactions suppressed by M and this leads 
to the inequality M 4 < A^ 5 //i 4 . This leaves room for interesting phenomenology as we will 



In (83) the factor of N e comes from the fact that the operator cr 4 keeps operating after 
modes have crossed the horizon. According to our prescription for the treatment of the con- 
version of cr fluctuations into density perturbations described in sec. 4, the non-Gaussianities 
induced when the modes are outside of the horizon are incorporated in the coefficients of 

15 For the UV oriented reader, it should be particularly easy to UV complete models with these kinds of 
symmetries, as the Lagrangian for the tt sector can be taken as the one of standard slow roll inflation. In 
general it is interesting to try to find particular UV completions of the models described by the Effective Field 
Theory of Inflation. But this is not necesary to have a consistent and predictive theory. 
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eq. (61), which are kept generic. Here however we explicitly include this effect in the estimate 
just to stress that, in the absence of cancellations from the terms coming from reheating, we 
expect the generated non-Gaussianity to be proportional to N e . 

For 1 < fx 4 N e M 4 /(H 4 A 4 U:S ) < (n s - l)N e M 4 /(H 2 A 2 U:S ), this four point function is the 
leading signal. It induces a qnl of the local kind of the order of 



a2 £4 
9nlQ ~ 7- 



Ne ~ ~^-N e , (84) 

E~H A U,S 



which can be detectable. A first analysis of this kind of four-point function in the WMAP 
data has been recently performed in [41]. In principle, the cubic operator in (28) of the 
form jl 4 ^ / A\j S is allowed by this symmetry and would give a larger effect on the three-point 
function. However, we saw that there are several symmetries that can forbid such an operator. 
On the other hand, the cubic operator in (26) schematically of the form a(d^a) 2 is negligible. 
By performing a similar comparison, we obtain 16 
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soft— breaking 
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soft— breaking 
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This could be only marginally detectable if we happen to see a local four-point function 
close to its current upper limit. We therefore conclude that we can have a detectable four- 
point function with the shape induced by the quartic operator a 4 . This shape is of the local 
form and is degenerate with what can be produced at reheating, but here we notice that 
this mechanism allows for having a large four-point function of the local form without at 
the same time having a local three-point function. In fact, as we just argued, we expect 
that non-linearities at reheating will induce non-Gaussianities dominated by the three-point 
function. 

There is another way in which we can obtain a large four-point function without having 
a detectable three-point function. We can impose the parity symmetry 

01 ->• -aj , (86) 

as an approximate symmetry in the Lagrangian. This symmetry is expected to be broken 
at the time of reheating. This symmetry removes all the cubic operators in a as well as the 
mixing ncr, and leaves only the cubic operators of the form 1x00 and the quartic of the form 
ixixua and <j<jo<j. Now, if the influence of n on ( happens to be much smaller than the one of 
the ct's (see eq. (72)) then the operators containing it are negligible and the leading interaction 
terms are of the form (da) 4 as contained in (20). Alternatively, we can impose the additional 

71 ->■ -7T , (87) 



16 Notice that in principle the soft-breaking operators like the one we consider here could appear in the form 
// sin(a-/A[/,s)(9(T) 2 /A^, as {da) 2 is compatible with the shift symmetry. For simplicity we consider only 
the case in which all terms are suppressed by A^ s . It is easy to check that our estimates would not change 
relevantly. The same applies to similar operators in the rest of the section. 



39 



approximate symmetry in the Lagrangian, and we are left again only with the (da) 4 operators. 
Because of the kind of symmetry we are imposing, the speed of sound of the a fields need not 
be equal to one nor equal to one another. For simplicity we will take all the sound speeds to 
be equal, the generalization being straightforward. In this case the speed of sound is given by 



l + e 2 



where we have suppressed the 1 index in because we are taking all of the speeds of sound 
equal. If all the coefficients Mi for the various operators are the same, in the presence of 
a speed of sound different from one different operators would contribute differently to the 
non-Gaussianities, their contribution differing by powers of l/c s according to the number of 
spatial derivatives versus time derivatives they have. However unless the coefficients of the 
various operators are constrained to be in some particular relationship by some symmetry it 
is reasonable to argue that there is a natural scaling relationship among the coefficients Mj 
such that they differ by powers of l/c s . The natural relationship is such that all of these 
operators become strongly coupled at the same energy scale, and therefore because they all 
have the same scaling dimensions, they have comparable effects at energy scales of order H 
that are relevant for inflationary perturbations. By natural relationship among the various 
coefficients we mean that if the coefficient of one of the operators does not respect this scaling 
relationship, then loop corrections will tend to renormalize with power law divergencies the 
various operators. This renormalization becomes an order one effect if the coefficients respect 
the natural scaling relationship. This statement is the exact analogous of the one we usually 
make for relativistic theories, when we say that all the various operators are suppressed by 
the same mass scale. In the presence of a speed of sound different from one, additional powers 
of c s need to be added. 

It is straightforward to find this natural scaling among the various coefficients. Since we 
are dealing only with the self-interactions involving the a-fields, there is no way for them to 
realize that their speed of sound is different from one: this means that we can redefine the 
spatial coordinates and the fields in such a way as to make the a fields have an effective speed 
of sound equal to one. Let us therefore redefine the spatial coordinates in the following way 

x -> x = — , (89) 

Cs 

and canonically normalize the a field as 

a ->. cr c = c\ /2 a . (90) 
In this way the quadratic Lagrangian becomes 

j 'dtd 3 x^^(& 2 -c 2 (<V) 2 ) -> Jdtd 3 x (a c 2 -(«V c ) 2 ) , (91) 

where <9j = d iv The effective speed of sound in this rescaled coordinates has become equal to 
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one. The quartic Lagrangian becomes 
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where in the first passage we have simply defined the M a ^ c in terms of the M\^,z of equa- 
tion (20). Since we did not rescale time, we can read off the unitarity bound directly 
by inspection of the mass scale suppressing the operators in the rescaled Lagrangian in 
(92). In particular, notice that if the M a fe c satisfy the following approximate relationship 

Afj ~ M c c 3 ~ M b c s ~ M a /c s , then all of these operators become strongly coupled at the same 
unitarity bound scale A^. It is arguable that loop corrections indeed force this relationship 
in a natural way. This line of reasoning was developed initially in [3, 5]. Unless explicitly 
specified in the rest of the paper we will assume this scaling in c s . As a result the g^L that 
can be generated by these quartic interactions is of order 
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where we have used (81). The induced gNL could be detectable. In this case there are three 
different operators that can induce a large four point function, 
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whose coefficients can be read off in (20) and are basically independent. Because these are 
derivative operators, the shape of the induced non-Gaussianity is of the equilateral kind, 
where all the sides of the quadrilateral in momentum space are comparable. Notice that only 
the shape generated by the operator a A is equal to the one that can be generated in single-field 
inflation [5]. The remaining three shapes can only be generated in multi-field inflation. 

The a A interaction coming from the soft symmetry-breaking terms can lead to a detectable 
local four-point function. It gives: 
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If we happen to detect this four point function close to its current experimental bound, it is 
possible that we will have the chance to detect also subleading operators, such as aa 3 , which 
is suppressed with respect to the leading term by a factor of H/A U:S and which gives rise 
to the same local four-point function. It is interesting to point out that in the inflationary 
case it is possible that symmetry is explicitly broken in a way that goes to zero in the 
limit the time-traslations are recovered as an exact symmetry 17 . In this case the lowest 

17 For example if the clock field were to be a standard scalar field, it is possible that the symmetry breaking 
is proportional to the inflaton time-derivative. 
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dimension operators associated to the soft symmetry breaking might be absent, subleading 
operators with derivatives acting on the cr's might become the leading ones, and operators 
that are not detectable if the lower dimension operators are present, such as those of the form 
fi A a 2 (da) 2 /Afj S , can indeed become detectable. Notice also that if the leading operator has 
only one derivative acting on the it's, then the upper limit on /i 4 that is determined by the 
induced mass of the cr's goes from /i 4 < (n s — l)if 2 A^ 5 to /i 4 < (n s — l)HAfj S . If instead 
the leading operator has two derivatives, then we have /i 4 <C Afj^ s . For simplicity in doing 
our estimates we shall consider the first upper bound /i 4 < (n s — 1)H 2 A 2 JS) but it is worth to 
keep in mind that these other possibilities are present. 

So far in this section we have concentrated on the case in which the dispersion relation 
of the a fields is of the form oo oc k. As we mentioned in sec. 3.1, the dispersion relation at 
horizon crossing can be of the form u 2 ~ k 4 /(e^M 2 ) if the parameter e 2 satisfies e 2 > M 2 / H 2 . 
This very non-relativistic dispersion relation is the same as the Ghost Condensate [38], and 
similarly to that case [39, 1, 3], we can have large non-Gaussianities even in this regime. The 
discussion proceeds in a very similar way to the case of a linear dispersion relations. While 
in this case we cannot impose that the theory be quasi Lorentz invariant in order to have a 
large four-point function without a detectable three-point function but we we can still impose 
an approximate a — > —o symmetry. Now the scaling dimensions of the operators are very 
non-relativistic. If we rescale t — > s~ 1 t, in order to keep the kinetic term invariant we need to 
rescale the spatial coordinates as x l — >■ s~ 1 ^ 2 x\ Invariance of the quadratic Lagrangian than 
implies that we need to rescale a as 



The leading quartic operators has scaling dimension 1/2 and is of the form (<9j<r) 4 . Remarkably, 
it is unique. This operator induces a four-point function of the order of 



for how this scaling is derived). Contrary to the case of single-clock inflation, here there is 
in general no symmetry protecting the generation of a standard spatial kinetic term (<9j<x) 2 
from a loop constructed by contracting two fluctuations of the (<9j a) 4 operator. Imposing 
that such a generation is still subleading with respect to the tree-level value of (di<j) 2 restricts 

~ = 1/8 ~ 

Mi < M/e 2 ■ Since the leading interactions are proportional to Mi, imposing such an 
inequality is radiatively stable. Notice however than since e<i > M 2 /H 2 , this constraint 
implies < H. Given the estimate in (97), this tells us that the theory with a (c^a) 4 
interaction has to be strongly coupled at energies scales of order H, and therefore not viable. 

There are two ways to proceed at this point. The first is to impose a symmetry that ensures 
that the operator (ftcr) 4 is small enough not to renormalize the (ftcr) 2 term. This can be done 
by remembering that as shown in [5], in single field inflation with the dispersion relation oj oc k 2 
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it is possible to impose an approximate n — > —it symmetry that makes quartic operators of the 
form 7r p (didjTr) 4 ~ p the leading interaction operators. This is done by keeping the contribution 
from the operators n 3 , ix(diix) 2 and (<9j7r) 4 small and by imposing their coefficients have a 
particular relationship among themselves so that their renormalization to the kinetic term 
(<9j7r) 2 cancel each other. This is the result of a symmetry. The renormalization of (<9j7r) 2 
needs to be proportional to terms that breaks the shift- symmetry of 7r and it is therefore 
small. For similar reasons the operators rr 3 , 7r(<9;7r) 2 and (<9j7r) 4 do not get renormalized by 
higher dimensional ones. If we impose that the a Lagrangian has the same structure as the 
7r Lagrangian then we can have the same signatures. In particular we can have a leading 
four-point function of the form & p (didjfj) 4 ~ p with size of order 
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where we have used that he unitarity bound Ay is given by A v ~ e5 /2 Mf /5 /M 3 / 5 . 

A second way to make the four-point function the leading signal in the case of uj oc k 2 
dispersion relation is by making the interactions associated to the soft-symmetry breaking 
the leading ones. In this case, the leading operator is the one in a 4 . Because of the non-linear 
dispersion relation, the constraint that the induced mass be smaller than ~ (n s — l) l l 2 H 
becomes 

H A < K - l)H 2 A 2 s . (99) 

Here we are normalizing the fields in the soft symmetry breaking sector by As, i.e. the 
soft-breaking potential takes the form ~ fi 4 cos(<t c /As). Because of the non-linear dispersion 
relation A5 is not equal to the unitarity bound Ajj t g in the soft-breaking sector but is related to 
it A| ~ A UtS M 3 e 3 2 /2 . As in the case of a linear dispersion relation, interaction terms compatible 
with the Abelian symmetry do not renormalize the soft breaking terms. However, the soft- 
breaking terms can generate interactions compatible with the symmetry. In particular the 
coefficient Mf suppressing the operator (did) 4 has to be smaller than Mf < A us M 5 e^ 2 //2 4 . 
They also generate the (pi a) 2 operator. Imposing that the generated term to be smaller than 
the tree-level value amounts to an inequality that it is satisfied for A^s ^ ( n s — l^^H, a 
constraint that we will easily be able to satisfy. The a 4 operator gives a local four-point 
function parametrically of the form 

£ (4) n 4 / TT \ 1/2 

soft-breaking ^ ^ |X < ^ _!)(») . (100) 



The operator da 3 , if present, would give rise to an additional contribution to the local four- 
point function suppressed by a factor of H/As and could be detectable. The operator a 3 , 
if present, would give a dominant and detectable local three-point function. The operator 
a(dicr) 2 is subleading with respect to the four-point function but is potentially detectable: 
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where in the next-to-last passage we have used that Ajj,s ^ e 2 M to avoid superluminality. 

It is easy to estimate that the quartic operator (<9j<r) 4 as generated by loops of the soft- 
symmetry terms gives a four-point function that is suppressed with respect to the one from 
o" 4 by a factor of H 2 /(Afj s N e ) and would be hard to detect. Instead, the operator a 2 (di<j) 2 
generates a four-point function that is suppressed with respect to one of a 4 by a factor of 
H I (A s N e ) that could allow for a subleading detection. 



• Non-Abelian case 

The story in the non-Abelian case proceeds in very similar terms for both the possible kind 
of dispersion relations. For simplicity we will mainly discuss the case of a linear dispersion 
relation. Given the Abelian treatment we have just studied the extension to the non-Abelian 
case is straightfroward apart for a small subtlety we will comment about. If we impose an 
approximate Lorentz invariance for the theory of the fluctuations, we remove all the cubic 
interactions and we are left with operators of the form 

en d ll a a d li a a d v a h d v a h + c 12 d^a a d'"a b d u a a d u a b . (102) 

Notice that at quartic order there are other less irrelevant operators (i.e. of dimension six 
instead of dimension eight) coming from the terms in D 2 in (38) that we wrote down in (44). 
If the mass scale F 1 suppressing these operators were the same as the one suppressing the 
operators in (102) then these operators would be much more important at energy scales of 
order Hubble, by a factor of order F 2 /H 2 ^> 1. However as we noted in eq. (65) through 
a rotation of the unbroken group, the adiabatic fluctuations are sensitive only to one linear 
combination of the a a fields. Without loss of generality it can be chosen to be the a = 1 
component. If are were looking only at adiabatic fluctuations then we would be looking at 
a four-point function of the <r a 's where all the indeces a are equal to one. Because of the 
antisymmetry of the structure constants C a bc, the contribution of the dimension six operators 
vanishes. For adiabatic fluctuations we are therefore left with the interactions in (102) that 
are of the same in form as in the Abelian case. 

As we noted in the paragraph below eq. (68) isocurvature fluctuations, if present, will in 
general be sensitive to a different linear combination of the <T a 's raising the possibility that 
the operators in (44) are important in the case of mixed adiabatic-isocurvature four-point 
functions. We will explore this further in the next subsections. 

The g NL induced by the operators in (102) is of the order of 
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and could be detectable. 

Operators similar to the ones in (102) are included also in terms of the form Tr [^D^^D^] 
and others with similar index contractions. Notice that naively these terms contain operators 
that are cubic in the a's and that are not removed by imposing a quasi Lorentz invariance 
for the theory. However, it is worth to check that they cancel because of the antisymmetry 
of the structure constants. 
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Another interesting case in which the four-point function is relevant observationally is 
when Ti[x a x a Xa\ = 0, which happens for some groups like a broken SU{2). In this case the 
three point function for the a fields is suppressed while the four-point function is not. This 
is another way in which a four-point function becomes the leading non-Gaussianity that it 
valid only in the non-Abelian case and that does not involve the requirement for the theory 
to be quasi Lorentz invariant. 

A detectable four-point function can also be associated with the soft-symmetry breaking. 
For example, the non-Abelian analogous of the operator in (23), schematically of the form 
a 4 , can give a detectable four-point function for \x in the interval 1CT 4 < /i 4 A r e /A 4 /5 < 
(n s — l)N e H 2 /Afj S . The non-Abelian analogous of the terms in (26), when present, could 
also be potentially detectable. The first term, schematically of the form a 2 &, can give a 
comparable effect in the three-point function. The second term, schematically of the form 
a(dcr) 2 , is suppressed by a factor H/(A UtS N e ), which could still leave room for a possible 
detection. It is interesting to point out that it is possible that the lowest dimension operators 
are absent because of the way the symmetry is broken. In such a case, subleading operators 
might become the leading ones, and operators that are not detectable if the lower dimension 
operators are present, such as those of the form n A a 2 (dcr) 2 / Afj s , can become detectable. 

Finally we add just one comment on the case the dispersion relation is of the form u oc k 2 . 
As we said, in this case everything proceeds as in the Abelian case except that the symmetries 
of the Lagrangian do not seem to allow to impose the same structure as in single field inflation 
and that the operators that generate the kinetic terms also induce some self-interactions. A 
straightforward estimate shows that these operators would renormalize the standard spatial 
kinetic term (di<j) 2 unless the unitarity bound is smaller than A v < F/e^ 6 . Here F plays 
in the non-Abelian case, the same role as M does for the Abelian one. While, as in the 
Abelian case, we cannot have interactions of the form (<9j<r) 4 , this inequality can be satisfied 
by introducing soft-breaking interactions whose unitarity bound A^g is smaller than F /e^ G . 
In other words, in the non-Abelian case with the dispersion relation uj oc k 2 , it is necessary 
to introduce soft-breaking interactions. This tells us that detection of a three-point function 
as produced by the operators a 4 , (dfcr) 4 and similar index contractions would rule out the 
possibility that the additional fields are the Goldstone bosons of a non-Abelian group. 

There are finally other possibilities in which the four-point function is smaller than the 
three-point function but still detectable that we will discuss in the next subsection. 



• Supersymmetric case 



In the case where the lightness of the a fields is partially protected by an approximate super- 
symmetry, there are two quartic interactions that can be relevant. The first is the AW* 2 
interaction. This leads to a four-point function of the order 
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(104) 



where we have inserted the factor of the number of e-foldings iV e coming from the fact that 
this interaction acts from horizon crossing to the end of inflation. Remembering that the 
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lightness of the a field is controlled by A <C 1, this four point function would be detectable 
in the interval 1CT 4 < \ 2 N e <C N e . We saw that such a shape for a four-point function is to 
a very good approximation the local one, which is produced also from non-linearities in the 
relationship between ( and a. As we will see in the next subsection, in the supersymmetric 
case this four-point function is naturally accompanied by a three-point function of the local 
form whose size is parametrically larger. The fact that we have not seen yet such a local three- 
point function constraints the four-point function in (104) to be practically undetectable. 

If we instead impose an approximate shift-symmetry on the imaginary part of a, then a 
four-point function of the local form induced by the operator A 2 Im(o") 4 can be generated and 
can be the leading signal. 

It is also possible to have a large four-point function induced by the derivative operator 
a^d^ad^a* /M 2 2 coming from the Kahler potential. By imposing an approximate E — > — E 
symmetry it is possible to suppress the three-point function, and the resulting g NL is of the 
order 
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The a-mass that this operator induces can be estimated to be of the order of H 2 /M St2 , where 
the loop is cutoff at energies of order H because of supersymmetry restoration. So, we should 
concentrate on the interval 1 <C M 2 2 /H 2 < 10 4 for detectability. However, we saw that 
supergravity corrections induce an operator of the form <7 4 H 2 /M 2 2 . The four point function 
induced by this operator is of the form 
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K2 M s,2 



where the factor of N e comes from the fact that this interaction keeps acting for all the time 
the modes are outside of the horizon. This is detectable in the interval 10~ 4 < H 2 N e /Ml 2 < 1. 
The size of this four-point function is a factor of N e larger then the one in (105) originating 
from the operator <t 2 (<9^<7) 2 and it has a local shape. It is therefore hard that the shape 
originating from <7 2 (<9 M <r) 2 can be detected unless we happen to detect a very large qnl of 
the local kind. We would be tempted to declare that such a detection would be a signature 
of supersymmetry as an approximate symmetry during inflation, the caveat being that in 
this case the lightness of the a fields is guaranteed only by a tuning of order n s — 1, and it 
is therefore delicate to make claims about the signatures of these models that are based on 
using naturalness at the percent level. 



We conclude this subsection by summarizing that it is possibile for multifield inflation 
to generate a large and detactable four-point function even in the absence of a detectable 
three-point function. Contrary to the case of single field inflation with a linear dispersion 
relation u> oc k, where only one shape is possible here there are five shapes generated by five 
independent operators. In some cases the signal is either peaked on equilateral quadrangular 
configurations (for the operators involving derivatives on each fluctuation) or in more general 
configurations (for the operators without derivatives on each fluctuation). We will study in 
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detail the shape of the four point function in a subsequent paper [36]. One unique combi- 
nation of these operators is Lorentz invariant. A detection of any of the shapes that are 
not produced by single field inflation would be a clear indication of multifield inflation, and 
depending on the particular shapes and a possible parallel detection of a three-point function 
we could learn if Lorentz symmetry or even, with some further luck, super 'symmetry are ap- 
proximate symmetries for the Lagrangian of the fluctuations during inflation. We could also 
have indication of the fact that the signals might be originating from soft-symmetry breaking 
terms. An analysis of the four-point function in the CMB data in search for the signal from 
the above operators is in progress [37]. 



5.2 Three-point function 

In the former subsection we concentrated on the four-point function as multifield inflation can 
predict an interesting structure for this signal that is completely new. This obviously does 
not excludes the possibility that in different regions of the parameters space we can have a 
large and detectable three-point function. Let us start with the Abelian case. 



• Abelian case 

The action in (11) has all possible interactions of the form 7r 3 , 7r 2 07, vraf and erf subject only 
to the constraint that a derivative must act on each fluctuations and the operator must be 
rotational invariant. Since the conversion of the n and o/'s fluctuations into ( in eq. (61) 
is local in real space, this means that each of these operators will give rise to a three-point 
function whose size will depend on the particular coefficients of (61) and on the particular 
operator chosen from (11), but whose shape will be equal either to the one of n 3 or to the one 
of 7r((9j7r) 2 . These are the two kind of shapes for the three-point function that are produced 
in single field inflation when the Goldstone boson is protected by an approximate continuos 
shift symmetry. Exactly the same shape can be produced in multifield inflation. For example 
in [31] the possibility of these large interactions had been highlighted by studying the case 
of multifield DBI inflation, though in that case the coefficients of a 3 and &(dicr) 2 are related 
because of the additional symmetry imposed on the fields by the DBI Lagrangian. In general 
in the Abelian case they are independent and they give rise to non-Gaussianities with size of 
order 

Jnl(, ~ 7- 

where now Ajj ~ cJ^M\ and can potentially be detected with high signal to noise. Here in this 
subsection we are assuming that the unitarity bound associated to the cubic term is smaller or 
equal to the one associated to the quartic interactions compatible with the symmetry, so that 
the leading signal is indeed the three-point function. The linear combination of the induced 
shapes interpolates continuously from shapes where the signal is dominated by equilateral 
triangles in Fourier space, whose amplitude is characterized by f^ lh , to shapes where the 
signal is dominated both by equilateral triangles and by flat triangles (with opposite signs), 
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whose amplitude is characterized by f^ L h ° s '- A detailed and optimal analysis of any possible 
linear combination of these two shapes was performed in [3] on the WMAP data finding no 
detection but constraining the amplitude of the signal of these two shapes to be [9, 3]: 

-214 < fPf- < 266 at 95% C.L. , -410 < /° r f og ' < 6 at 95% C.L. , (108) 

with f^ L hog ' being zero almost excluded at 2a level. Using the effective field theory of single 
field inflation these constraints were then translated into contour plots for the coefficients 
M 4 and M| of the single field operators 7r 3 and 7r(<9j7r) 2 [3]. Those contour plots could be 
extended to include all the additional operators we have in the case of multifield inflation, 
though, given the large degeneracies (many operators give in fact only two shapes) a large 
fraction of the parameter space would be left unconstrained. We leave this for future work. 
Here we just mention that among the many different ways in which the three-point function 
can be large there are the cases in which the 07 fields have a small speed of sound and there is 
also the case in which the time-kinetic mixing is characterized by having the e unm i x defined in 
(13) much smaller than one. Of these two cases, the second one is an intrinsically multifield 
effect. 

Operators that induce a shape for the three-point function that is different from the ones 
produceable in single-field inflation are the ones associated to the soft-breaking of the U(l) 
symmetry of eq.s (25) and (26). First lets consider the case when the symmetry breaking is 
happens due to only one spurion, the leading operators are the following. The operator in 
(26), schematically of the form //cr 2 a- c /A| ~ ^Hal/ A A S) gives a local three-point function of 
size 
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where we have used that fi A < (n s — l)H 2 Au S , which could be detectable. This is the leading 
signal. Subleading signals are given by the operators fi 8 a 3 /A 7 s , whose signal is suppressed 
by an additional factor {n s — 1)H/Au,s, and /x 4 cr(9cr) 2 /A|, whose signal is suppressed by an 
additional factor H/(A U:S N e ). These signals could be marginally detectable. Further, as we 
saw in sec. (5.1), we can impose an approximate Lorentz invariance symmetry which forbids 
the cubic operators compatible with the U(l) symmetry as well as the operator schematically 
of the form a 2 a from symmetry breaking. In this case the leading cubic operator has the 
form a(d fl a) 2 . This is unique and it induces a shape for the three-point function that cannot 
be produced as the leading ones in single field inflation, and whose size is of the order 
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which can be detectable in the regime 10~ 4 < if/Z/A 4 ^ < (n s — l)H 3 /Afj S . Notice that the 
non-Gaussianities induced by this operator are expected to be smaller than the ones induced 
by the operator generating the local four-point function in (23), schematically of the form cr 4 , 
by a factor of H/ (A Ut sN e ). Because of this, detectability of such a shape is possible only in the 
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case we happen to detect a large four-point function of the local kind. Notice that while this 
shape could be produced in single-field inflation as well, it is in general either subleading with 
respect to different shapes than the local four-point function and/or accompanied by other 
shapes of comparable size. Therefore the combination local four-point function and shape 
of the form <j(<9 M o") 2 is a clear indication of multifield inflation with an approximate Lorentz 
invariance. Further notice that for the interactions of the form cr(dcr) 2 one of the fluctuations 
does not have a derivative acting upon it. Therefore the resulting shape of the three-point 
function will have a non-vanishing squeezed limit and so it will affect the clustering statistics 
of collapsed objects such as galaxies and there would be constraints in studies similar to [42] 
in addition to direct searches for a three-point function. 

In this same regime where we suppress the leading cubic operators by imposing an ap- 
proximate Lorentz invariance for the fluctuations, the local three point function induced by 
the soft-symmetry braking <r 3 operator in eq. (25) could be detectable. Its non-Gaussianity 
is suppressed by a factor of order < (n s — l)H/A UjS , which could still allow for a detection. 

In the case the explicit symmetry breaking is performed with two spurions of different 
charge such that a large cubic operator of the form <r 3 is generated, the induced local three- 
point function is of the size 
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which can be detectable in the interval 10~ 4 < N e j2 4 / (Ay S H) < if/A^g . We remark that it 
is possible that the soft symmetry breaking might be proportional to the parameter breaking 
time diffs. (related to the inflationary clock field). In this case, the leading soft-breaking 
terms might be absent, making the subleading ones even more detectable. 

As in the former section when we discussed about the four-point function in the Abelian 
case, we have so far concentrated on the case where the dispersion relation is linear in k, of 
the form u = c s k. As we discussed, it is also possible that the dispersion relation at horizon 
crossing is of the form u = k 2 /(e^ 2 M). Because of the non-relativistic scaling, in this case 
the leading operators are of the form a(dia) 2 and (d 2 a)(di<j) 2 . The level of non-Gaussianity 
induced by the first operator is of the form 
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where A^ ~ e^J 2 M s / M 7 . A similar estimate holds also for the operator (d 2 a)(di<j) 2 , and 
the resulting signal can be detectable. The same signal can be produced also in single field 
inflation, where optimal limits on the three-point function as induced from these two operators 
were given in [3] by deriving and interpreting the limits on f^ lh and /™2 h ° g '- 

In analogy to what happened for the four-point function in the absence of a symmetry the 
operator in a(di<j) 2 can generate at loop level a kinetic term of the form (di<j) 2 . Imposing that 
such a generation is still subleading with respect to the tree- level value of (<9j<r) 2 restricts M < 
M/e 2 • Since the leading interactions are proportional to M imposing such an inequality is 
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radiatively stable. The constraint e 2 > M / H 2 forces the unitarity bound A v to be smaller 
than H. Given the estimate in (112) such a situation is not viable. This would naively 
force us either to forbid the cubic interactions compatible with the symmetry, or to impose a 
symmetry by which we can make sure that the operator (di<j) 2 does not get renormalized. We 
can go back to the case of single-clock inflation in the near de Sitter limit with a dispersion 
relation of the form oo oc k 2 where there are three interactions of the form 7r 3 , 7r(<9j7r) 2 and 
(<9j7r) 4 with some specific coefficients fixed by symmetries. There can also be an additional 
unrelated interactions proportional to (<9 2 7r)(<9j7r) 2 , accompanied by symmetry reasons by 
higher order interactions. Contrary to the naive estimates loops of these operators cancel 
each other and do not generate a (diir) 2 term. This is in fact a result of a symmetry: the 
(<9j7r) 2 term has to be proportional to HM^ and such a term breaks the shift symmetry of ir. 
It can therefore be generated only through operators that explicitly break this symmetry. This 
ensures that there is a cancellation among the loop-diagrams which might seem remarkable 
at a diagrammatic level. Coming back to the case of multifield inflation where the a fields 
have a dispersion relation of the form w oc fc 2 , we see that if we impose that the coefficients 
of the operators ct 3 , a(di(r) 2 and (<9;cr) 4 to have exactly the same relationship as in the case 
of single-clock inflation, then we are guaranteed not to generate the (di<j) 2 term (and that 
the relative coefficients between the operator remains unchanged). In this way we are able 
to eliminate the constraint M < M/e^ 8 and we take M ~ M/e^ 4 so that the unitarity 

-1/2 = 

bound Ajj ~ e 2 M coincides with the scale at which a excitations become superluminal. 
The resulting non-Gaussian signal is dominated by a three-point function as in (112) with a 
subleading, and undetectable, four-point function as in (97). 

Let us now consider the symmetry breaking terms. In the case we break the symmetry 
with only one spurion, the leading operators are of the form /i 4 cr c <T 2 /A| ~ /i 4 if of/A 4 , and 
/i 4 <r c ((9jO" c ) 2 /A|. The first one gives rise to a local three-point function of the form 
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where we have used that Me^ 2 > A UyS , that A| ~ A Ut sM 3 e^ 2 and that /i 4 < (n s — l)if 2 A|. 
The second operator gives a three-point function of the order of 
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Notice that the three-point function induced by these two operators is smaller than the local 
four-point function induced by the a 4 operator, which scales as (n s — l)N e (H / A^s) 1 / 2 ! . Still 
the signal could be detectable even in the presence of the a 4 operator. 

If the symmetry is broken with a second spurion jl 4 , then it is possible to have a large 
local three-point function as induced by the /i 4 cr 3 /A| operator. The signal scales as be 
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which can be detectable and which, if present, would make the contribution from the other 
operators negligible. In order for the operators associated to soft-breaking not to renormal- 
ize the (did) 2 kinetic term, we have to impose Au,s ^ ( n s — l^^H, a constraint that is 
straightforward to satisfy. 

In multifield inflation, there is one additional form of three-point function that does not 
have an analogue in single field inflation. This is due to the three-point function induced by 
the quadratic relationship in (61) between ( and a. This induces a three-point function of 
the local form, with as in (74). The optimal analysis of the WMAP data for this kind 
of signal was performed in [9, 43], where the constraint were combined with analogous limits 
from the SDSS data [42] to obtain 

-5 < < 59 at 95% C.L. . (116) 

This limit can be translated, using (61), into limits on the parameters that characterize the 
conversion from 07's to (. Unfortunately, the large degeneracy in the parameter space limits 
the power of this constraint. 

It is worth stressing that single field inflation cannot produce a value of f 1 ^ that is larger 
than the deviation of the power spectrum from scale invariance, which is generically very 
small [17, 18, 2]. This means that a detection of such a three-point function in the absence 
of a detection of a large deviation from scale invariance would rule out single field inflation. 

• Non-Abelian case 

The non-Abelian case is a straightforward generalization of the Abelian one and for simplicity 
we shall treat explicitly only the case of a linear dispersion relation, adding just a brief 
comment concering the non-linear dispersion relation towards the end. The main difference 
is that the possible operators, and therefore the possible shapes, can be further constrained 
depending on the non-Abelian group considered. If Tr[:r a ;r{,:r c ] 7^ 0, then we can have the 
same kind of cubic interactions a 3 and &(di<j) 2 present in the Abelian and in the single-clock 
cases. Notice that because of the field redefinition (65), in the case of adiabatic fluctuations 
we are actually interested in groups where Tr[x a x a x a ] 7^ 0. We have cubic interactions of the 
form 

a 3 , &(8 t a) 2 , only if Tr{x a x a x a ] ^ . (117) 

The trace is non-zero whenever there is a commuting Abelian subgroup (which had to be the 
case as this case includes the Abelian case of the former section). But for other groups such 
as SU (2) it can be zero. 

When the trace is zero we are left with the last three cubic operators in (43) and with 
the terms from soft-symmetry breaking. Notice however that the last three operators in 
(43) require that both the o fields and the 7r field be important at a comparable level for 
the adiabatic fluctuations. The operators from soft-symmetry breaking require a relevant 
breaking at the time of horizon crossing, which again does not need to be there. We find that 
this strengthen the relevance of the four-point function in the non-Abelian case. 



51 



We now come back to the operators associated with the soft-symmetry breaking. The 
natural suppression of the three-point function for some non-Abelian groups makes the dis- 
cussion even more interesting than in the Abelian case. If we break the symmetry with only 
one spurion, the leading operator might be given by the non-Abelian analogous of the cubic 
one in eq. (26), schematically of the form a(da) 2 . This gives an effect of the order 
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In the limit 10~ 4 < /// (HAfj s ) <C H 3 /Afj s this signal is detectable. A simple estimate shows 
that the non-Gaussianity induced by this operator is subleading by a factor H/ (A UjS N e ) <C 1 
with respect to the four-point function induced by the non-Abelian analogous of the operator 
in (23), schematically of the form a A . There is however a regime in which both shapes 
could be detected. Notice that if the nature of the symmetry breaking is such that the 
first of the operators in (26) is present, then there is a local three-point function whose 
size is comparable to the one of the four-point function. The stronger constraints currently 
available from the data for the local three-point function [3, 9] imply that detection of the 
signal from (118) would be problematic. However, as we stressed, the presence of any of 
these operators depends on the characteristics of the symmetry breaking terms and it can be 
absent. This is for example the case for the well-known chiral SU(2) x SU(2) symmetry of 
the QCD Lagrangian spontaneously broken to the diagonal subgroup and explicitly broken by 
the quark masses. Here there are even more possibilities, as the explicit symmetry breaking 
could be proportional to the breaking of time diffs. by the inflationary clock, possibly forcing 
the leading symmetry breaking terms to have derivatives acting on the cr's. 

Because of the natural suppression of the three-point function that can happen for some 
non-Abelian groups the local three point function induced by the soft-symmetry breaking 
a 3 operator could be detectable. If the operator in aa 2 is present its non-Gaussianity is 
suppressed by a factor of order < (n s — 1)H/ A^s with respect to the one from the a 4 operator. 
This makes it comparable to the non-Gaussianity from the operator a(da) 2 in (118) and it 
could still allow for a detection. As in the Abelian case this operator can be made much more 
important and actually the leading one from the symmetry breaking sector if we break the 
symmetry with the non-Abelian analogous of two spurions with different charge. 

We finally add a brief comment on the non-Abelian case when the dispersion relation is 
of the form u oc k 2 . We notice that since the non-Abelian analogous of the kinetic operator 
in (14) induces automatically higher order terms, imposing the Lagrangian to have the same 
structure as in single-clock inflation to prevent the generation of the {di<j) 2 term seems hard. 
Non-renormalization of the standard kinetic term forces the unitarity bound to be smaller than 
F/e^ 6 , which is smaller than the scale ~ e^ 2 F at which these interactions become strongly 
coupled. As before we have to forbid the cubic interactions compatible with the symmetry, 
and introduce soft-breaking terms with a unitarity bound Au,s ^ F ■ Again soft-breaking 
interactions are in this case necessary. This also means that detecting a three-point function 
coming from the operators &(di<j) 2 or d 2 a(di<j) 2 with a dispersion relation ou oc k 2 excludes 
the possibility that the additional fields are Goldstone bosons of a non-Abelian group. 
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• Supersymmetric case 



In the supersymmetric case we can have two kind of three-point functions. The first comes 
from the operator Xm s (a + a*) a a* ~ XH(a + a*)aa*, where we have used that m s ~ H 
in order to cancel the mass from supergravity corrections. This operator induces a possibly 
detectable three-point function of the order of 



JnlQ ~ 7- 

J~2 



N e ~\N e =s> f NL ~ 10 5 AiV e , (119) 



which in the interval 10~ 4 < XN e < N e is detectable. The shape is local and the factor of N e 
comes as usual from the fact that this interaction keeps operating after horizon crossing. 

Another operator that can induce a large three-point function is the higher derivative one 
(a + a*)d ll ad^a* /M Sjl . This induces a three-point function of order 



InlC ~ 7^ 

J~2 



" h ^ w Wr (120) 



which would be detectable in the interval 1 <C M Sjl /H < 10 4 . The lower limit on M s l comes 
from imposing that the cx-mass induced at loop level is much smaller than one (after cutting 
off the loop at scales of order H where supersymmetry is restored). As we saw, such a term 
in the Kahler potential induces also a cubic operator of the form if<7 3 /M Sil from supergravity 
corrections. The induced three-point function is of the local form and has a size of order 



InlC ~ 7^ 



N *~l^- N e f NL ~10 5 ^-N e . (121) 

E~H M *,l 



This is a factor of N e larger than the one in (120) and is detectable in the interval 1CT 4 < 
HN e /M s> i < 1. Given the current constraints on /]^£ al , the signal in (120) will be hard to 
detect. 



5.3 Isocurvature fluctuations and their non- Gauss ianity 

As we anticipated in sec. 4.1 in the case of multifield inflation isocurvature fluctuations can 
be generated. In this case a fluctuations would be visible not only through their effect on the 
adiabatic fluctuations but also through their imprint in the isocurvature ones. Of course this 
has been known for a long time (see for example [47] and references therein for the two-point 
function, see for example [24] for non-Gaussian correlation functions). Our results build upon 
the former ones. We believe the framework we give in (68), together with the Lagrangian 
for the fluctuations, offers a general parametrization for the study of these effects and for 
the identification of new signatures. Furthermore it offers the most general way of relating 
experimental constraints (or measurements) to parameters of the theory. 

From the linear components of (68) it is straightforward to see that we can have isocur- 
vature two-point functions of various kinds, as well as cross correlations between any two 
kinds of isocurvature fluctuations and between adiabatic and isocurvature fluctuations. Cross 



53 



correlations happen in cases when two different fluctuations depend on the same fields or if 
they depend of two different fields that are mixed at quadratic level in the Lagrangian. This 
mixing can only happen between the n field and the Goldstone boson of a commuting Abelian 
subgroup (see sec. 3.1). 

Isocurvature fluctuations can be non-Gaussian. As in the adiabatic case there are sev- 
eral mechanism that can lead to non-Gaussianity. The simplest example is of course if the 
quadratic terms in (68) are relevant. This would induce an isocurvature three-point function 
of the local type. For the same reason as in the adiabatic case we do not expect a four-point 
function to be generically induced in this way. 

Alternatively, through the linear dependence in the conversion parameterization of (68), 
non-Gaussianities in the statistics of the a fluctuations can be translated into the isocurvature 
sector. Here the discussion develops as in the case of adiabatic fluctuations: we can have three- 
point functions with the same shape as in the adiabatic case. Furthermore by imposing the 
same kind of symmetries as in the adiabatic case to the multifield Lagrangian we can have a 
large four-point function without a detectable three-point function in an isocurvature mode. 
However, there are important differences with respect to the purely adiabatic case that we 
are now going to explain. 



• Abelian case 

For the Abelian case, the possible shapes for the three and four-point function arise from 
the same kind of operators (dcrj) 3 and {daj) 4 where there is a derivative acting on each 
fluctuation, plus the operators of the form cr(<9 M <r) 2 , <r 4 , bo 2 and a 3 , if present, associated 
with soft-breaking of the U(l) symmetry. Further, in the case of isocurvature fluctuations 
we can have more than one kind of three and four-point functions depending on the kind of 
correlation being considered. All the correlations of the form 



(CCCisoC. 



(122) 



are possible, where we have defined C iso as one of the isocurvature flucutations in eq. (68) 
and where the number of insertions of ( and Ciso is in general arbitrary. Since we know 
observationally that isocurvature fluctuations are smaller then the adiabatic ones, the leading 
effect comes from inserting the smallest possible number of isocurvature fluctuations in the 
correlation functions. The [/(Incompatible operators give rise to a non-Gaussianity of the 
order of 



Co 



E~H 
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(123) 



Notice that when the four-point function is made the leading signal by imposing some symme- 
try of the a Lagrangian, the scale suppressing the cubic and the quartic operators is different, 
the scale associated to the quartic operators being actually lower in order for them to be 
relevant (see former subsections). Here we have for simplicity used the same scale. Here we 
have also defined e iso as the suppression of the isocurvature fluctuations with respect to the 
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adiabatic ones and we have inserted only one factor of e iso . WMAP data already constrain 
the isocurvature component to be subdominant so e iso < 0.1 — 0.3, the actual limit depending 
on the kind of isocurvature fluctuations considered [9, 47]. Forecasts for Planck give roughly 
e iso ~ 10~ 2 [46]. The ?7(l)-soft-breaking terms of the form a 2 &, a{d ll a) 2 and a 4 , whenever 
present, instead give a level on non-Gaussianity of the form respectively 
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(124) 
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where the factor of N e comes from the fact that the a 4 and do 2 interactions keep operating 
outside of the horizon. The a 3 operator, if present, would give a level of non-Gaussianity of 
order 

£3 



~4 -4 

TV ~ f . — — — TV ~ f- — — — TV 

Jy e t iso. 3 rr jy e t iso. 3 rr Jv e> 



(125) 



which would be the leading one. 

1 /2 

The limits on the ratios in for example (123) and (124) scale as 1/iV 4 , where Ap ix is the 
number of signal dominated modes. Clearly if ei so is not too small non-Gaussianities involving 
isocurvature fluctuations might be detected. 

As in the former subsections, so far we have dealt with the case in which there is a linear 
dispersion relation for the modes. In the case the dispersion relation is of the form oj oc k 2 , 
the discussion proceeds in very analogous terms. The analogous of Eq. (123) becomes 
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with the value of Ajj given in the former subsection. 

The terms associated to the explicit symmetry breaking in eq. (124): a 2 a, a(di(r) 2 and a 4 
give: 
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where in the last passage we have used that e 2 ' M > A^g. Given the current constraints 
on the non-Gaussianities, all of these terms could be detectable. If the term in a 3 were to 
be present with a spurionic coefficient proportional to /i 4 , it would induce a level of non- 
Gaussianity of order 
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which could be detectable as a local three-point function. 



• Non-Abelian and Supersymmetric cases 

Everything in the former discussion generalizes to the non-Abelian case and the supersym- 
metric case. There are however some important differences concerning the three and the 
four-point functions in the non-Abelian case that we now discuss. For the three-point func- 
tion we noticed in the adiabatic case that there are additional cases where the three-point 
function was naturally suppressed when only one field was relevant for the adiabatic fluc- 
tuations. This was because we required that Tr[x a x a x a ] ^ 0. Observation of a three-point 
function in this case required that the explicit symmetry breaking was large enough. In the 
case of isocurvature fluctuations it is somewhat more natural to expect that more than one 
field is relevant for observations, and that might alleviate the suppression of the three-point 
function. 

A more interesting interaction is the one in (50) of the form a(da) 2 that can appear in the 
non-Abelian case when there is explicit symmetry breaking and one of the Goldstone bosons 
gets a time dependent vev. In this case the level of non-Gaussianity for a three-point function 



is approximately of the order of 



18 
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~ e iso — — ~ e iso — . (130) 
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This combination of isocurvature and adiabatic fluctuations is chosen in order to minimize 
the suppression from inserting isocurvature fluctuations. This signal is expected to be dom- 
inated by at least an adiabatic four-point function of the local form, that is larger by a 
factor of iV e /ei SO , and so could be detectable only if we happen to detect a large four-point 
function. If the operator associated to soft breaking a 2 a is also present, it induces a local 
three-point function in the adiabatic sector comparable to the four-point function. Given the 
more stringent constraints we have on the local three-point function, in this case detection 
of the signal in (130) would become problematic. As we stressed, the presence of the various 
symmetry-breaking operators in the non-Abelian sector depends on the spurionic transfor- 
mation properties of the terms explicitly breaking the symmetry. It might be that some 
operators are absent in specific cases. 

Concerning the four-point function, we saw in eq. (44) in sec. (3.2) that there are dimension 
six operators of the form a 2 (da) 2 which could not contribute to the four-point function of 
adiabatic fluctuations because of the antisymmetry of the structure constants and to the 
fact that adiabatic fluctuations could depend only on one linear combination of the <r a 's 
(see eq. (63)). For the adiabatic case this forced the leading operator contributing to the 
four-point function to be of dimension eight. As we pointed out in the paragraph below 
eq. (68), the linear combination of fields that sources isocurvature fluctuations is in general 
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See footnote (16) for an explanation on how the fields are normalized in this case. 
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different from the one sourcing the adiabatic fluctuations. This implies that the operators 
in eq. (44), schematically of the form a 2 (da) 2 , can give a non-zero effect in the case of 
isocurvature fluctuations. Of course in this case the effect is suppressed because we know that 
observationally isocurvature fluctuations are smaller than the adiabatic ones. For this reason 
the leading contribution will come from mixed adiabatic-isocurvature four-point functions of 
the form 

(CCCisoCiso) • (131) 

which is non-zero and has the smallest suppression from the smallness of the isocurvature 
fluctuations. The amount of non-Gaussianity induced by the operators in (44) schematically 
of the form a 2 (da) 2 is of order 
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Notice that in the adiabatic case, the level of non-Gaussianity scales as 
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The ratio of these two four-point functions scales as (e iso A;y/iJ) 2 which, if isocurvature fluc- 
tuations are not too small, might lead to the leading detectable signal. Notice that these 
dimension six operators are absent in the Abelian case. Detection of such a kind of mixed 
four-point function with such a shape in the absence of a large non-Gaussianity in the adi- 
abatic fluctuations is therefore possible and it would be a clear indication of a non-Abelian 
sector of Goldstone bosons active during inflation. 

The same shape can also appear in a correlation function of the form 

(CCCCiso). (134) 

This is possible when there is explicit symmetry breaking and the non-Abelian Goldstone 
bosons develops a time-dependent vev. In this case the operator in (51) again of the form 
a 2 (da) 2 gives a non-Gaussianity of order 



u 

^iso r 
1-2 



~ CisoTT-Ts fiso-T^ , (135) 



which, if // happens to be close to its upper bound (n s — l)if 2 Aj^ 5 and e iso not too small, 
could be detectable. This same shape as induced by an operator of the form a 2 (da) 2 can 
also be generated in the supersymmetric case. However in this case we would not expect a 
suppression of the non-Gaussianities in the adiabatic sector and, as we have seen, we would 
expect an even larger three-point function of the local form. 

Another interesting shape proper of the non-Abelian case that can appear in ((( CCiso) is 
given by the operators in (45), schematically of the form a(da) 3 . They induce a four-point 
function that scales as 

e iso -w- ~ £ iso-T^ ~ e isoTy • (136) 

^2 Er ^ H Ajy 
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Since one of the fluctuations does not have a derivative this shape will have a non- vanishing 
squeezed limit and might be detectable by studying the clustering statistic of collapsed objects. 
Because of the non-Abelian symmetry, these interactions have to be dominated by cubic ones 
of the form (da) 3 by a factor of Fi/(e iso H). In turn they give a level of non-Gaussianity that 
scales as 

If we detect a large level of non-Gaussianity in the three-point function and if e; so is not 
too small then we might marginally detect this four- point function. Such a correlation of 
a signal between a three-point and a four-point function would be a clear indication of the 
non-Abelian nature of the symmetry group. 

The situation when the dispersion relation is 00 oc k 2 is also very interesting. In the three- 
point function the effect of the operator considered in (130) schematically of the form a(dia) 2 
becomes 



^iso r 
>~>1 



H 4 ( H \ 1/4 
~ e iso = - , /s < e iso (n s - 1) . (138) 

Though in the non-Abelian case one might not have all the naively allowed symmetry breaking 
terms, it is worth noticing that if an adiabatic interaction of the form a 4 is present, then the 
induced signal would compare to the above one by a factor of N e (H / Au^s) 1 ^ (with possibly an 
additional factor of 1 / ej so in case the a 4 operator contributes in the purely adiabatic sector), 
which could be either dominant or subleading depending on parameters. The same signal 
would be suppressed by a factor of N e (modulo an additional possible factors of l/e iso as 
before) if the operator of the form a 3 were to be present with coefficient jl 4 . 

For the operators considered in (132) schematically of the form a 2 (da) 2 , we notice that 
the leading term is the one with two time derivatives a 2 a 2 . While naively, because of the 
non-linear dispersion relation of the form u oc k 2 , the term in a 2 (dia) 2 would be much more 
relevant (it would actually be a dangerous relevant operator), the structure of the Lagrangian 
forces its coefficient to be suppressed by a factor of l/e 2 < (H/Au) 2 with respect to the 
coefficient of a 2 a 2 which therefore becomes more relevant. Its effect scales as 
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with F playing the same role in the non-Abelian case as M in the Abelian case and we have 
used that F > A U)S e^ 6 . 

Considering now the operator in (135), schematically of the form a 2 (d i a) 2 ) its effect scales 

as 

u 



^iso ~ 
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<e iso (n s -l) — (140) 



e~h iS ° HV 2 Af s M^~eT ~ " S A%t 



This compares to the local four-point function by the factor M 1 / 4 e^ 8 /(A^A r e ) which could 
be either larger or smaller than one. Detection of such a shape would be another stricking 
signature of the non-Abelian symmetry group. 
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The spatial-kinetic operator in the non-Abelian case, schematically of the form (T>iD l a) 2 , 
induces interactions of the form a 2 (dfa) 2 where one of the fluctuations has to be isocurvature. 
The induced four-point function is of order 
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which is subleading with respect to the leading ones. 

Finally the operator considered in (136), schematically of the form cr(da) 3 , would give a 
four-point function of the form 



This is subleading with respect to the three-point function that has to be present at the same 
time by a factor of order €\ so (H / hu) 1 ^ , which would make it very hard to detect. 



Given the ongoing experimental effort to test inflation and the proliferation of different models, 
it is quite important to characterize the most general theory of inflation. In this paper we 
have extended the recently developed effective theory of single-clock inflation [1, 2, 3] to 
include additional light scalar fields that are in their vacuum at horizon exit. In a subsequent 
paper [12] we shall include non-scalar light fields and fields that are not in their vacuum state 
at horizon exit. By concentrating on the theory of the fluctuations, we have been able to 
write down the most general Lagrangian coupling the Goldstone boson of time-translation (71) 
with additional light scalar fields. We have assumed the additional scalar fields to be either 
the Goldstone boson associated with the spontaneous breaking of some symmetry, or to be 
protected by an approximate super symmetry. These two mechanisms provide a natural way 
to make scalar fields light. 

Contrary to the case of single-clock inflation, the Lagrangian for the fluctuations around 
the time during which the observable modes cross the horizon is not sufficient to make full 
contact with observations. In principle, knowledge of the full Lagrangian in field space is 
necessary to determine how fluctuations in the additional light fields get converted into per- 
turbations of cosmological interest. However, most of this conversion happens when gradients 
are negligible and this allowed us to parametrize the general conversion mechanism with a 
few parameters, without significantly affecting the predictive power of the effective theory of 
multifield inflation. 

We have been able to identify many observational signatures that can be produced both 
in single-clock inflation and in multifield inflation, such as a large three-point function of 
the density fluctuations. However, we have found a very rich structure of possible additional 
signatures specific to multifield inflation, such as some specific shapes of the three-point 
function and more interestingly, the realization that it is quite generic for multifield inflation 




(142) 



6 Conclusions 
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to produce a large and detectable four-point function with a larger signal-to-noise ratio than 
the three-point function. This opens up a new observational window into inflation that has so 
far been largely unexplored. We have also shown that there are specific signatures associated 
to the Abelian or the non-Abelian nature of the symmetry group of which the additional light 
scalar fields might be the Goldstone bosons off, as well as signatures associated to a possible 
approximate supersymmetry of the Lagrangian during inflation. Depending on the specific 
signatures we might detect, we could be able to distinguish between these three mechanisms 
for keeping the scalar fields light. Our estimates for the size and the characteristics of the 
signals we have identified have been approximate and valid at the order of magnitude level; 
clearly a more detailed study on the lines of the one performed in [3] for single-clock inflation 
should be performed in the future. 

Concluding, we can summarize the two main results of our paper by saying that we have 
developed a universal framework where all multifield inflation models can be studied at the 
same time and their signatures analyzed, and we have realized the importance of a new, 
largely unexplored possible signature of inflation: the four-point function. A detection of 
such a signature would teach us a lot about the dynamics during inflation and even possibly 
about the fundamental symmetries of Nature. 
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